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CHAPTER  I 

;  .  ■  INTRODUCTION  ari  rar  -^ 

An  important  problem  in  highway  design  is  that  of  determining  flow    '  ' 
capacities  required  for  drainage  structures .    Although  the  problem  exists 
for  structures  of  all  sizes,  there  is  a  partic\ilar  shortage  of  design  criteri 
for  small  structures  such  as  culverts «    Culvert  installations  ordinarily 
are  used  where  the  discharge  originates  from  small  watersheds  of  a  few 
acres  or  a  few  square  miles »    The  determination  of  peak  discharge  magnitudes 
and  corresponding  return  frequency  intervals  is  essential  to  economical 
engineering  design o 

Until  the  year  1955  there  was  very  little  reliable  information  available 
regarding  -peak  discharge  magnitudes  and  frequencies  from  small  drainage  areas 
in  Montanac    In  that  year  the  UoSo  Geological  Survey  instituted  a  Small- 
Area  Peak-Flow  Hi^way  program  in  cooperation  with  the  Montana  State  Highway 
Commission o    The  cooperative  program  involves  the  determination  of  peak 
flow  rates  at  selected  stations  on  small  watersheds o    Measurements  were 
begun  at  ^5  crest  sta^e  gaging,  .stations  in  1955,  with  more  stations  being 
added  in  1959  and  later  years,  with  228  stations  being  in  operation  in  1969« 

In  order  to  derive  the  maximum  possible  utility  and  benefit  from  the 
cooperative  progrmn  a  comprehensive  study  of  peak  flows  from  small  watersheds 
(l  to  100  square  miles)  in  Montana  was  undertaken  by  Montana  State  University 
in  1963o    Principal  objectives  of  the  study  were  to  obtain  information  which 
would  at  least  partially  answer  the  following  questions:  -  »■ 


lo     Can  existing  precipitation  and  climatological  data  be  lised  to 
predict  the  frequency  of  flood  magnitudes  on  small  watersheds  in  Montana. 


2o    Do  caaprehensiTe  studies  of  a  few  of  the  watersheds  in  the 
peak-flow  stu*^  of  the  UoSo  Geological  Survey  (ioeo^  detailed  precipitation  and 
snow 'measurements 9"  soil  types g  topography ^  land  usej  continuous  flow  measure- 
ments   etco )  yield"  data  which  ejctend  the  usefulness  of  the  results  obtained 
from  the  peak  flow  determinations  being  made  by  the  Geological  SurTeyo 

3o    'What  frequency  of  flood  can  be  safely  predicted  from  short-term  peak 
flow  records  0    What  length  ©f  records  should  be  made  by  the  Geological 
Surrey  to  safely  predict  floods  of  given  frequencies  0 

J".  ?■! 

In  the  furtherance  ©f  the  Project  objectives ,  the  study  was  initially 
divided  into  two  phases i 

1  -  Precipitation  data  gathered  by  the  UoSo  Weather  Bureau  (ESSA) 
were  assembled  and  studied  to  determine  what  correlation  exists  between 
such  data  and  peak  flows  from  small  watersheds  o    In  this  phase  data 
already  available  or  being  obtained  by  others  were  utilizedo 

2  CoH^rehensive  hydrologic  studies  of  four  of  the  UoSo  Geological  Survey 
watersheds  were  madeo    Data  obtained  were  correlated  with  the  results  of  the 
Geological  Survey  programo    In  this  phase  new  data  were  gathered  specifically 
for  this  Projeeto 

As  the  investigation  progressed^  the  scope  of  the  Project  was  enlarged 
to  include  a  study  of  various  methods  current^  in  use  by  other  agencies  for  the 

-  2  - 


prediction  of  peak  flawB  and  frequencies ,  ond  this  led  to  a  critical  assess- 
ment of  all  modem  hydrologic  teciiniq.iies  and  their  applicability  to  Montanao 

The  varioms  phases  of  the  investigation  are  reported  in  separate  chapters 
of  this  report o    Reference  is  frequently  made  to  the  one  preliminary  report 
and  seven  interim  reports  which  have  been  prepared  and  submitted  during  the 
life  of  the  Project o    Although  it  is  intended  that  this  final  report  Tbe  complete 
and  self-contained a  many  details  which  appeared  in  the  earlier  reports  are 
only  "briefly  smrmariaed  herein o 

,  A  review  of  modem  hydrologic  techniques  appears  in  Appendix  Ao       •  -•: 

Several  graduate  students  and  a  large  mmber  of  undergraduate  students 
at  Montana  State  University  were  employed  at  various  times  on  this  Project o 
The  conscientious  efforts  of  these  individuals  should  be  acknowledgedo  Faculty 
colleagues  s,  Dro  Eo  Ro  Dodge  of  the  Department  of  Civil  Engineering  and  Engineering 
Mechanics^  Dro  To  Lo  Hanson  of  the  Department  of  Agricultural  Engineering, and 
Dro  Ko  Jo  Tiahrt  of  the  Department  of  Mathematics  served  as  technical  consultants 
on  the  Projects,  and  eaceh  contributed  immeasurably  to  the  study o  Professor 
Ao  Co  Seheer  of  the  Department  of  Civil  Engineering  and  Engineering  Mechanics 
was  Project  Coordinator g  rmd  advised  on  fiscal  and  administrative  matters  o 
Dro  Paul  Brown 9  Agricultural  Research  Services,  USDA^  obtained  the  loan  of  an 
Infiltrometero    Technical  personnel  fi^om  the  UoSo  Geological  Survey ^  ESSAg 
and  Soil  Conservation  Service  were  consulted  freely,  and  without  exception 
they  were  most  cooperative o    A  total  of  l8  ranchers  in  eastern  Montana  not 
only  permitted  the  inetallai'fcion  of  meteorological    instruments  on  their  property, 
but  serviced  the  equipment  regularly  and  mailed  in  charts  and  records  faithfully 
each  montho    Finally  inention  should  be  made  of  the  Montana  Highway  Commission 


and  UoSo  Federal  Highway  Administration  personnel j  both  in  Helena  and  at  District, 

offiees  throughout  the  state  j  who  aided  the  Project  in  comitless  respects., 


( 


CHAPTER  II 

v.-    ^  ■         PRECIPITATION-STREAMFLOW  CORRELATIONS 

Phase  I  of  the  Proposal  for  this  Project  envisioned  the  compilation  and 
study  of  precipitation  data  gathered  by  the  UoS»  Weather  Bureau  (now  the 
Environmental  Science  Services  Administration)  to  determine  what  correlation 
exists  between  these  data  and  -peak  flows  on  small  watersheds «    This  chapter 
describes  the  procedures  that  have  been  followed  in  pursuance  of  this  objective. 

ESSA  operates  nearly  350  precipitation  measurement  stations  in  Montana 
of  which  about  90  are  equipped  with  continuous-recording  weighing  raingages , 
More  than  100  of  the  stations  have  been  in  operation  at  least  50  years o  No 
significant  amount  of  precipitation  data,  except  for  that  at  the  ESSA  stations, 

was  discoveredo  -..v  '       '  ^  .--a' 

Interim  Report  #7  (Robinson  and  Williams,  1968)  describes  the  development 
of  a  method  to  relate  peak  discharge  to  maximum  rainfall  intensity  on  small 
Montana  watersheds o    Although  data  available  are  not  sufficient  to  adequately 
test  the  method,  results  obtained  are  encouragingo    It  is  believed  real  ■ 
possibilities  exist  for  development  of  more  rational  culvert  design  procedures c 

RAINFALL  FREQUENCY  -  PEAK  FLOW  FREQUENCY  METHOD 

In  relating  frequency  of  rainfall  intensity  to  frequency  of  peak  flow 
it  is  first  recognized  that  a  peak  flow  which  has  a  certain  return  period,.?;-, 
(say  ten  years)  may  be  produced  by  any  one  of  n.umerous  rainstorms  Depending 
upon  a  host  of  antecedent  moisture  factors,  watershed  characteristics  and 
meteorological  conditions,    the  ten  year  discharge  might  be  produced  by  a 


six-hour  diaration  rainstorm  of  one-inch  per  hour  intensity,  by  a  three-hour 
duration  rainstom  of  two  inches  per  hour  intensity,  or  by  a  twelYe-hour 
duration  rainstorm  of  one-half  inch  per  hour  intensity »    If  soil  conditions 
were  rights  the  same  ten  year  discharge  could,  in  fact  be  produced  by  a 
six-hour  duration  rainstorm  of  one-half  inch  per  hour  intensity,  or  it 
might  require  a  six-hour  duration  rainstorm  of  two  inches  per  houro    In  other 
words  an  infinite  variety  of  rainstorms  might 5  given  the  proper  conditions, 
produce  a  peak  discharge  of  ten  years  return  periodo 

>nr-  J.:.. 
Although  there  would  seem  to  be  little  possibility  of  a  simple 
function  relating  10  year  peak  flow  to  10  year  rainfall  intensity  (ioeo- 

-  K  X  I^^q)  it  could  be  hypothesized  that  at  higher  return  periods  such 
a  relationship  would  become  more  valido    This  is,  a  rainfall  of  such 
intensity  that  it  occurs  once  in  50  years  would  be  esqpected  to  produce  a  peak 
discharge  having  a  return  period  close  to  50  years;  and  even  more  extreme 
events  might  be  expeeted  to  follow  even  more  closely  the  flmction       -  K  x  I^ 

The  prediction  equation  which  was  developed  is  based  upon  the  above 
hypothesis  and  is  given  by 

^  Q  ^  (1) 

where       is  the  peak  discharge  of  return  period  i 

is  a  rainfall  intensity  ratio j  being  the  ratio  of  the  i-year 
rainfall  intensity  of  some  given  duration  to  the  mean  annual 
rainfall  intensity  of  the  same  duration 


~  6  - 


is  a  rainfall~dis charge  reciirrenee  factor  which  expresses  the 
relationship  between  the  rainfall  recurrence  relation  for  specified 
duration  and  the  recurrence  relation  of  peak  annual  ^schargeso 
is  a  rain-snow  baseflow  interaction  ratio  which  reflects  the  relative 
dependence  or  independence  of  the  frequency  curves  of  rainfall- 
induced  flows  and  snowmelt-induced  flows  o  ■■  ;\ 
Q  is  the  mean  anniaal  peak  discharge  rate  of  the  streanio 

The  return  period  i  was  taken  as  50  years  in  this  study c 

Deterzaination  of  R. 


Precipitation  data  from  l8  ESSA  stations  in  eastern  Montana  with  recording 

raingages  and  26  stations  with  non-reeording  raingages  were  examinedo  Peak 
2^-hour  intensity  was  foimd  for  each  year  for  the  non-recording  stations, 
and  peak  one-houTj  peak  two-hour,  peak  three-hour 9  peak  four-hour »  peak 
six-hour^  peak  eight-hours,  peak  12-hour ^  and  peak  2^-hour  intensities  were 
obtained  for  each  year  for  the  recording  stations o    These  data  were  analyzed 
statistically  (by  Gimbel^s  technique)  to  determine  intensities  corresponding 
to  50-year  and  to  2o 33-year  return  periodso    The  ratio  of  these  two  values,   :  ,  ■ 
I^q/I  was  determined  in  each  case 9  and  called  the  rainfall  intensity  ratio. 

It  was  found  that  for  rainstorm  duration  greater  than  12  hours  5  the  •, 
values  of  R  were  remarkably  consistent  (see  Interim  Report  #7?  Figure  11, 
page  66)0    This  led  to  the  conclusion  that  R  values  could  be  obtained  from 
non-recording  station  data  md  these  values  could  be  applied  to  watersheds 
having  travel  times  in  excess  of  12  hours c 


A  map  of  Montana  showing  isoplethals  of  R^^  was  shown  in  Interim  Report 
#7  (Figure  12^  page  68)0    Since  preparation  of  the  report  further  work  has 
been  done  on  refinement  of  R  values ^  and  the  map  which  is  shown  herein  as 
Figure  1  is  believed  to  more  aceurately  define  R^^  in  eastern  Montanao 

Determination  of  D„  .  . 


Yaltaes  of  the  rainfall-diseharge  recurrence  factor       were  obtained 
for  eleven  watersheds  in  eastern  Montana  varying  in  area  from  30  to  6Qk 
square  miles ^  having  continuous  streamflow  records  varying  from  11  to  36 
years o    In  analyzing  streamflow  records  from  these  watersheds  only  that  portion 
of  the  streamflow  assumed  attributable  to  rainfall  was  usedo  Snowmelt- 
ihduced  runoff  was  separated  out  of  the  hydrographs  where  necessary  by  stu<fying 
characteristic  shapes  of  those  hydrographs  which  evidently  were  totally 
caused  by  snowmeltj  and  subtracting  comparable  flows  from  hydrographs  which 
were  partly  rain"=caused  and  partly    snowmelt-inducedo    Baseflow  was  separated 
out  of  the  hydrographs  by  assming  that  the  baseflow  equalled  the  average 
of  the  recorded  mi.nimum  daily  discharges  during  the  month  in  question □ 
was  computed  for  each  watershed  from  eactreme  value  distribution  theory, 
and  is  given  by 


1  C 


1 


(in  i  »  Y^^) 


(2) 


1  +  c 


1 


where  C    is  coefficient  of  variation  (standard  deviation  divided  by  the 


V 


mean)  for  discharge  peaks  (Q)  or  maximum  rainfall  intensities  (l). 


(p    BSk3    y>,  are  escpected  values  of  standard  deviation  and  mean  for  the 
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length  of  record  of  discliarges  (Q)  or  intensities  (l)o 
i    vm  taken       50  yeaxso 

Isoplethals  of  D^^  for  the  area  of  Montana  east  of  the  continental  divide 

were  shown  in  Interim  Repor'fc  #7  (Figure  13 ^  page  T0)o 


Determination  of  F^^ 


ValTues  of  the  rain-snov  baseflwdatermeti to:  ratio,  F^^^,  were  computed 
for  eleven  watersheds  (those  -ysed  in  determination  of    D^g,  described  above) 
from  the  relation 

;     ^1  =  ^L/A.  (3) 

where  Qj^  and  Q^^  are  pesQt  diseharge  at  i=-year  return  period  and  mean 

annual  discharge j  due  to  rain  only 

'  i 

and  Q^^  are  peak  discharge  at  i-year  return  period  and  i 
meaK  annual  dis charge o 
i  was  taken  as  50  years o 


Values  of  F^^  were  found  to  be  quite  erratic,  and  not  apparently 
related  to  geographical  or  climat ©logical  factors o  F^^  was  therefore 
tentatively  taken  m  lo^5  (the  average  of  the  values  obtained  for  the  eleven 
watersheds ) o 
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TEST  OF  METHOD 


A  50 -year  peak  discharge  was  eomputed  by  the  rainfall  frequency-pealt 
flow  freqiaency  method  for  eaeh  of  k6  of  the  watersheds  equipped  with  crest- 
stage  gages o    These  peak  diseharges  were  eompared  to  eomparable  values  determined 
hj  the  extreme  Talue  funetioa  (Gmbel  method)  and  by  the  log-probability 
aiethodo    Results  ©f  the  eomparisons  are  shown  in  Figures  2  and  3°    For  most 
©f  the  k6  watersheds  the  three  methods  predict  very  nearly  the  same  dis charge o 
Qe  a  few  watersheds  there  is  eonsiderable  discrepancy^  In  most  of  these 
oases  the  method  developed  herein  gives  smaller  diseharges  than  either  of 
the  methods  lased  for  cheeking o  :  •  ■ 

DISCUSSION  OF  RESULTS  -• 

There  iSg  of  course ^  no  final  means  of  evaluating  the  rainfall  frequency- 
peak  flow  fret^ency  method  (or  any  other  method  for  that  matter)  because  no 
values  are  available  for  cosnaparisono    This  method  has  one  important  advantage 
over  other  methods ^  in  that  it  is  less  sensitive  to  length  of  streamflow 
reeordo    Stresimflow  records  are  used  only  to  obtsdn  mean  annual  discharge o 
The  variance  from  the  mean  is  obtained  from  precipitation  records  which 
are  usually  much  longer  than  the  streamflow  records c    A  disadvantage  is  that 
each  ©f  the  three  parameters       D  and  F  must  be  developed  separately  for  each 
recurrence  frequency  desiredo    R^^  for  instmce  is  different  from  R^.  and 
one  esyriaot  be  eoi^uted  from  the  other o  .  ■ 

laterim  Report  #7  (Robinson  and  Williams »  1968)  describing  the  rainfall 

fr@queney=peafc  flow  frequency  method  was  reviewed  by  Do  Co  ¥oOg  for  the  Structural 

and  Applied  Meehaaies  Division^  Federal  Highw^  Admini strati on o    The  following 

coraaents  are  intended  to  partially  answer  some  of  the  questions  raised  by  WoOo 
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FIGUEE  2.     Comparison  of  estimates  of         between  the  peak 

rainfall  intensity- peak  flow  method  and  extreme 
value  frequency  analysis. 


-  12  - 


1830 

N) 

—      r  ■■    ■  :  •  .y- 

1209 

1771  5 
O 

0  o 

\  ^ 

1832 

20.317 
\^  cfs 

1765 

\ 

1395 

W  1 

1 

1 

2  4  6  8 

Estimated  by  Log- Probability  Analysis  (cfsxlO  ) 

50 

FIGURE  3.     Comparison  of  estimates  of         between  the  peak 

rainfall  intensity- peak  flow  method  and  log- 
probability  analysis.  ,^  ,.. . 
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lo    Woo  "Wondered  whether  the  extreme  value  distribution  (Gumbel's 
technique)  Qm  properly  "be  applied  to  rainfall  intensities  in  Montanao 
(Gumbel's  teehnique  was  used  to  determine  precipitation  intensities  having 
a  JO-year  return  period^  and  hence  to  determine  values  of       =  R^q/Ro) 
Confidence  eurves  have  been  drawn  for  the  2i|-hour  precipitation  intensities  at 
13  of  the  stations  used  in  the  study  and  also  for  the  one-hour,  two-hour, 
three-hours,  four-houTs,  six-hour,  eight-hour,  12-hour  and  2i+-hour  precipitation 
intensities  at  Martinsdaleo    These  curves  have  been  plotted  on  the  Gumbel 
frequency  graphs  for  these  stations  and  are  shown  in  Figures  k  to  2ko  The 
confidence  curves    show  the  range  of  return  periods  within  which  a  given 
magnitude  of  the  precipitation  intensity  may  be  expected  with  a  probability 
of  about  2/3o 

I 

Figures  k  to  2k  also  show  90  percent  confidence  bounds  on  the  mean 
annual  precipitation  intensity,  using  the  assumption  that  the  population  of 
rainfall  intensities  is  a  normal  distribution,  and  calculated  by  Student's  t, 

i 

(x  -  /a)  ^nr 

s 

where  t  is  Student's  t 

(two-»tailed  test,  90%  probability  used) 

X  s  saaple  mean 

population  mean 

S  ~  sample  standard  deviation 

n  ~  number  years  record 
Also  shown  on  Figures  h  to  2h  are  90  percent  confidence  upper  bounds  on 
the  $0-year  precipitation  intensity,  calculated  by  again  assuming  that 
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preeipitatioE  intenslt^r  lAieh  is  within  the  2/3  confidence  limits  'utsing 
Gmbel^s  teehnig.ii@9  hut  lower  thaa  Gijmbel"s  50-year  Talueo 


Study  of  FigiDires  k  t©  2k  leads  to  tl 
teotoit^e  may  Tbe  siaeeess fully  applied  to 


Moat  ana o 


Wo©  '^uek^jtioaad  whether  the  rainfall  fret^eney-peaJi  flow  frequency 


¥SB  inde©d  pri 


SIS  aire 


The 


Qm  we 5 


Kiel 


¥alu@s  of  F  are  all  much  lower  than  the  "sigaifis^t"  the 


IS  tna 
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The  F  gtatistle  is  obtained  as  the  ratio  of  treatment  mean  square  to 
error  mean  sqiuareo 

Treatmeat  meaa  s^mare  ^  treatment  swa.  of  squares /degrees  of  freedomo 
Error  mean  sqiiaare  ^  error  bwi  of  squares /degree  of  freedom 


Tre@,tment  degrees  of  freedom  =  t  ~  1 
Error  degrees  ©f  freedom  =  t(r=-l) 

3c„  =  a  value  of  50=-year  peaJi  discharge  at  station  i  computed  "by  one 

of  the  3  methods  o 

r    "  niifflilber  of  itations  (kS  in  this  ease) 
t    -  amber  of  methods  (either  2  or  3  in  this  ease)o 
3o    ¥©©  noted  that  the  method  cannot  be  verified  because  there  are 
no  eacact  eorreet  values  to  be  xmed  for  comparison o    There  can  be  no 
argment  with  Wo©°s  statement  5>  but  it  must  also  be  applied  to  the  Gnambel 
teehniqiie^  the  lOrg^-normal  method^  the"  Log"°Pearsoa  Type  III  method,  or 
indeed  to  my  method  ^lich  purports  to  predict  the  magnitude  of  a  flood 
peak  having  a  ispeeified  return  peri  ©do    The  contention  that  the  method  is 
not  usable  because  of  the  absence  of  actual  long  records  is  not  believed 
to  be  valido 

¥©0  pointed  ©'at  that  the  values  of       D  and  F  used  for  estimating 
peak-discharges  were  derived  from  drainage  areas  of  30  to  TOO  square  miles,  and 
that  the  method  should  n©t  be  used  to  predict  peak  flows  from  drainage 
areas  outside  this  remge  of  areas  c    This  is  a  valid  point  j  and  should  be 
adhered  tOo 
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TABLE  I.     Comparisons  of  Estimates  of         by  Three  Methods 
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125.7 
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Wets  Creek  8.1 
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Pryor  Creek         |  435 
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12 
12 
12 
42 
19 


105.7 
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83.3 
82.0 
64.8 


1240 
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369 
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932 


1090 
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2570 
859 


Extreme  Value  (Gumbel)  Estimate. 
''^Log- Normal  (log  probability)  Estimate.  (Chow) 
'^Pec-k  Rainfall  Frequency- Peak  Flow  Estimate. 
^Adjusted  for  the  Event  of  June  1964. 
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TABLE  I.     Comparisons  of  Estimates  of         by  Three  Methods.  (Cont'd) 


8 


O 


U.S.G.S. 

Station  No.  • 

Watershed 

Area,  sq.  mi. 

Years  Record 

Mean  Annual  F 
(cfs) 

Coefficient  < 
Variation  {^0 

Estimate  c 
Q 

cu 

)f  50-year 
(cfs) 

0 

rd 

Peak 
0 

a 
0 

Ui 
•H 

0 

3082  1  Basin  Creek  Trib 

3083  j  Basin  Creek 
32^7      Sand  Creek 
3329    I  North  Creek 
33^1    '  Wolf  Creek 

0.1 

10.9 
10.6 
0.7 
9.1 

12 
12 
12 

13 
12 

52.2 
4l2 
56.2 

312 
396 

209.8 
92.6 
134.8 

104.0 
75.9 

439 
1756 

324 
j  i403 
1456 

351 
1545 

258 
1283 
1274 

205 
1620 

224 
1310 
1610 

3364.5  Spring  Creek 

3.9 

11 

111 

102.8 

I  515 

454 

590 

Extreme  Value  (Gumbel)  Estimate. 
''^Log-Normal  (log  probability)  Estimate. 
Peak  Rainfall  Frequency- Peak  Flow  Estimate- 


^  X 

fex)% 


87908.8 

352593000 
167999000 


74901.1 

618457000 
121960000 


59580.8 

225599000 
77170800 


D.F. 

SUM  SQUARES 

mm  SQUARE 

F 

^(.01) 

Gumbel- 
C  how- 
Robinson 

Among  Methods 
Within  Methods 

2 

135 
137 

8,742,000 
1,329,520,000 
1,338,262,000 

4,371,000 

9,848,290 

.443834 

4.61 

Gumb  el- 
Chow 

Among  Methods 
Within  Methods 

1    1  1,839,130 

90  11,181,090,000 

91  11,182,930,000 

1,839,130 
13,123,200 

. 140143 

6.97 

Gumb el- 
Robinson 

Among  Methods 
Within  Methods 

1    i  8,722,620 

90  i  833,022,000 

91  !  841,745,000 

1 

8,722,620 

9,255,800 

.942395 

6.97 

Chow- 

Among  Methods 

1 

2,551,250 

2,551,250 

.356029 

6.97 

Robinson 

Within  Methods 

90 

644,925,000 

7,165,830 

91 

647,476,000 

—  -  _ 
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CONCLUSIONS 

The  rainfall  freq'y.ene3r=>peak.  fl©¥  frequenej  method  is  an  important  contri- 
bution to  small  watershed  hydrologjo    For  a  watershed  with  at  least  a  short- 
tem  streamflow  reeord  available  it  cgyn  be  applied  rapidly  by  simple  formulas  <> 
It  does  not  reqmre  e3ctensive  field  datao    It  is  unique  among  methods  reported 
in  the  literature  in  that  it  utilizes  long-term  precipitation  data,  which  are 
available  in  qumtityj  to  determine  varianee  fr-m.  'mm  annual  flood  peako  It 
has  been  shpwn  to  give  results  that  are  eomparable  t©  those  obtained  from  other 
methods o 

Further  work  is  needed  to  refine  the  proeedure  and  eactend  it  to  other 
geographieal  areas  md  other  return  periods  o 

The  development  of  a  rainfall  frequency-peak  flow  frequency  method  for 
predicting  peA  discharges  ©f  given  recurrence  interval  is  discussed  herein » 
The  statistical  comparisons  and  tests  which  have  been  made  are  believed  to 
saaswer  the  principal  objections  raised  by  the  Federal  Highway  Administration 
reviewer  of  Interim  Report  #7o 


.    •    .    CHAPTER  III 
HIDROLOGIC  STUDY  OF  FOUR  SELECTED  WATERSHEDS 


Phase  II  of  the  Proposal  for  this  Project  envisi  oned.  the  compre- 
hensive hydrologic  study  of  four  widely  separated  watersheds  for  the 
purpose  of  extending  the  usefulness  of  the  U^So  Geological  Survey  "Small- 
Area  Peak-Flow-Highway  Progranio"    This  chapter  defines  the  scope  of  the 
studies  that  have  been  madeg  describes  the  procedures  used,  and  discusses 
data  collection  and  reduction  techniques  and  problems o 

SELECTION  OF  WATERSHEDS  ■ 

The  first  criterion  for  selecting  the  watersheds  that  were  to  be 

studied  was  that  they  must  be  watersheds  that  had  already  been  included 
in  the  UoSoGoSo  investigations  and  for  which  at  least  a  few  years  of 
crest-stage  gage  measurements  were  available o    In  conference  with 
Hydraulics  Division  personnel  of  the  State  Highway  Department,  and 
representatives  of  the  Uo So  Go  So  and  UoSo  Weather  Bureau,  some  guidelines 
for  watershed  selection  were  set  forth o    The  area  west  of  the  continental 
divide  was  eliminated  because  relatively  few  flood  problems  on  small 
watersheds  west  of  the  divide  could  be  citedc    It  was  agreed  that  all  the 
watersheds  should  be  in  plains  environments,  and  several  general  areas  of 
the  state  were  siiggestedo 

Thirteen  watersheds  were  visited  and  inspected  before  the  four  study 
watersheds  were  seleetedc    Brief  descriptions  of  the  watersheds  visited  were 
included  in  Interim  Report  #1  (Williams,  1965)0    Maps  and  aerial  photographs 
of  the  watersheds  selected  are  presented  in  Appendix  B 


Bacon  Creeks  in  Wheatland  County  drains  21o2  square  miles,  and  is  a 
tributary  of  Antelope  Creek  which  in  turn  is  trilbutary  to  the  Musselshell 
Riverc    Bacon  Creek  watershed  is  an  almost  exeltisively  range  area  (mostly 
sheep  grazing)  due  north  of  Harlowtono    Drainage  is  predominantly  north- 
west to  southeast'o    There  is  essentially  no  timber  on  the  watershedo  It 
has  had  a  history  of  occasional  extreme  flood  events j  the  largest  event  of 
record  having  "been  in  Juneg  1950  o 

Duck  Creek ^  in  Prairie  and  McCone  Counties  drains  ^hoQ  square  miles, 
and  is  a  tributary  of  Redwater  River o    The  watershed  has  two  main  streams. 
East  Duck  and  West  Duckj,  which  merge  a  short  distance  above  the  mouth  of  ' 
the  watershedo    Drainage  is  northerly o    The  watershed  extends  south  to 
Sheep  Mountain 3,  which  is  the  divide  between  the  Yellowstone  and  Missouri 
drainages  0    The  higher  elevations  are  rough  and  broken j  similar  in  character 
to  mmy  Badlands  areas  o    This  portion  of  the  watershed  is  used  almost 
exclusively  for  range  and  livestock  grazing o    Near  the  mouth  of  the  water- 
shed slopes  are  gentle  to  rolling  and  a  considerable  portion  of  the 
area  is  in  dryland  wheat o 

Evmp  Creek „  in  Sweetgrass  Coimty,  drains  7° 6  square  miles,  and  is  a 
direct  tributary  of  the  Yellowstone o    Drainage  is  to  the  northo    Mouth  of  the 
watershed  is  two  miles  west  of  Reed  Point  o    This  is  the  smallest  of  the 
watersheds  seleetedg  and  is  the  only  one  with  significant  growths  of 
timber 5    (mostly  scrub  pine)o    The  basin  is  steep ^  and  is  characterized  by 
numerous  side=-canyons  which  are  tributaiy  to  the  main  streamo    The  basin 
is  mostly  covered  with  native  grass  and  is  used  for  livestock  grazingo 

Lone  Man  Coulee,,  in  Pondera  County,  drains  lkc6  square  miles,  and 

is  located  five  miles  south  of  Ifalier^  and  about  two  miles  south  of  Lake 
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Frances  c    The  watershed  is  tributary  to  the  Dry  Fork  of  the  Marias  River, 
Drainage  is  from  west  to  easto    Slopes  are  gentle  and  most  of  the  basin 
is  in  cultivated  strip-^wheat  fanningo    The  largest  flood  of  record  at     •  - 
Lone  Mmi  Coulee  vae  a  flow  of  1820  cfs  in  June  19^8 >  hut  this  was  almost 
matched  by  a  flow  of  IT^^O  cfs  in  June  196h.    The  196h  flood  was  the  subject 
of  ^  speeial^  .interim  report  o  (Williains,  196^)0  y^-A'-.  , ,y.'^ 

PHYSICAL  CHAMCTERISTICS  ■       ■-     ■  -  ■  ' 

A  coiH.pils.tion  of  physieal  characteristics  of  the  project  watersheds 
is.  in  Appendix    G      <,    Included  are  data  related  to  area,  shape,  elevation, 
slope stream  leng'fch,  stream  density,  orientation,  pondage,  and  land  use. 

aerial  photographs  and  field  studies  were  utilized 

in  the  ciet@nd,nation  of  the  parameters .    The  soil  investigations  in 
particular  vere  qmte  esc&ensive  and  will  be  discussed  separately. 

Soils  liwestigaition  •  ■■■■■■■  ■      '  :•■  .  ■ 

Soil  analyses  on  the  four  project  watersheds  included  mechanical  analysis 
and  Atterberg  Limit  determinations;    medium  intensity  soil  sxirveys  by 
soil  seientiBts  of  the  Soil  Conservation  Service,  USDA;  and  water  intake 
studies.  ■   V  ■  .  '■■   .  ■ 

Meehanical  Analysis  and  Atterberg  Limits:    Soil  samples  were  obtained 
from  S'-ineh,  9°ineh  md  l8-inch  depths  at  each  of  the  eight  weather  station 
instnsEent  enclosures  at  the  time  of  installation  of  the  instruments  in 
196^0    The  depths  were  selected  to  coincide  with  the  depths  specified  for  soil 
moisture  and  soil  temperature  element  installation o  •  ■  ■ 
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Combined  sieve  and  hydrometer  analysis  were  made  on  each  sample,  sieve 
sizes  included  1-inch 5  ih^  #20  and  0hO^  with  the  hydrometer  being  used  on 
the  -^0  fractioBo    Particle  diameters  as  small  as  OoOOl  millimeter  were 
d@t@rHiined«.-: ,    b^j:--'  .'x;:'   :  -    .^■ 

Liquid  limit  and  plastic  limit  were  determined  for  each  sample,  and 

plasticity  index  was  calculatedo  ^  - 

Bmrnary  soil  data  are  shown  in  Appendix  D 

Soil  Surveys;    The  Montana  State  Office  of  the  Soil  Conservation  Service 
provided  soil  maps  and  descriptions  of  each  of  the  four  watersheds,  under  a 
cooperative  financial  agreement  with  the  project o    None  of  the  project 

watersheds  had  been  surveyed  for  soil  elassifieation  previously o  About 

90  days  of  field  work  were  required  by  the  Soil  Scientists  on  the 

four  watersheds  after  which  the  SCS  Cartographic  Unit  in  Portland,  Oregon 

prepared  the  maps  0    The  state  soil  scientist  personally  toured  each  watershed 

with  the  mapping  scientists g  md  gave  close  supervision  to  the  preparation 

of  final  maps  and  reports o 

-.The  SCS  classes  the  surveys  made  on  the  watersheds  as  "medium  intensity 

soil  surveys  in  which  strongly  contrasting  soils  are  separated  and  where 
areas  of  15  acres  or  larger  are  mappedo"    Soil  types  delineated  on  the 
finished  maps  range  in  number  from  12  each  at  Hump  Creek  and  Lone  Man  Coulee 

to  3^  at  Duck  Creeko  The  maps  are  all  prepared  on  a  scale  of  IS2O5OOO  (3<.l68 
inches  »  1  mile) o  ■ '  j  ■ 


Tabulated  in  Appendix  D      are  summary  sheets  for  each  watershed,  showing 
soil  types  delineated 3  approximate  extent  of  the  watershed,  and  SCS  hydrologic 

classificationo  ......  ... 

Water-Intake  Studies,     Infiltration  capacities  were  determined  in  I966 
for  the  major  soil  types  found  on  the  project  watersheds,  using  a  "mobile 
raindrop  applicator"  which  was  loaned  to  the  project  by  the  Agricultural  • 
Research  Service,  USDAo    This  device,  which  is  pictured  in  figure  25  ,  supplies 
simulated  rainfall  at  a  uniform  rate  to  a  circular  area  of  approximately  13 
square  feet  from  a  height  of  ih  feet  above  the  groundo    Water  runoff  measure- 
ments are  made  from  a  two-foot  square  test  plot  located  in  the  center  of 
the  area  receiving  rainfallo    With  this  arrangement  the  test  plot  is  relatively 
free  from  the  influence  of  lateral  water  movement o 

■  ■  ; 

The  raindrop  applicator  was  developed  by  Mto  Frank  Rauzi,  Soil  Scientist 
for  the  Agricultural  Research  Service  at  Laramie,  Wyoming^    It  was  made 
available  for  this  project  through  his  generosity,  and  the  efforts  of  Dr. 

Paul  Browng  Soil  Scientist,  ARSj,  Bozemmic    Raus^i  has  published  several 
articles  describing  his  results  with  the  raindrop  applicator  on  Northern  ^ 
Plains  watersheds o    His  procedures  and  methods  were  discussed  in  a  196O 

article 0  •■■^••^  -  ■ 

Test  runs  of  60  minute  duration  were  made  at  each  of  69  locations  on 
the  four  watersheds =    The  Soil  Conservation  Service  soil  scientist  assigned 
to  each  watershed  was  consulted  as  to  the  most  representative  locations  of 
the  various  soil  types  and  infiltration  tests  were  made  at  the  locations 
suggested  by  this  scientist  o    At  le^t  two  tests  were  made  at  each      ,   ■  ■ 
site,  to  check  reproducibility o    Tests  were  also  made  near  each  of  the 
eight  weather  station  installations o 


Mobile  Raindrop 


FIGURE  25  —  Mobile  Raindrop  Applicator. 


-  k6  - 


A  test  consisted  of  applying  rain  at  the  rate  of  three  inches  per  hourg 
with  nm-=off  from  the  test  plot  being  recorded  at  five  minute  intervals  o  Water 
intake  w©^  measured  as  the  difference  between  applied  rainfall  and  measured 
runoff  o  ■ '         '  '      '  '  .^.  r 

Soil  samples  for  moisture  determination  of  the  0-6  and  6-12  inch  ■ 
soil  depth  were  taken  near  each  test  plot  iaimediately  prior  toj  immediately 
after 9  and  2^  hours  after  each  testo    Percentage  soil  moisture  was  determined 
for  each  sample  o  ■•    -  ■  '  ■ 

All  standing  vegetation  and  mulch  material  in  the  test  plots  were 
clipped  at  ground  level  after  the  testg,  and  air-dry  weights  determinedo 

Depth  of  penetration  of  water  was  determned  at  each  test  plot  ■■ 
immediately  after  a  test  ^  and  again  2^  hours  afterwardo  •        .■  ,  ■ 

Water  int^e  rates  were  plotted       function  of  time  since  start  of 
testo    Simmary  sheets  for  each  watershed  are  in  Appendix     D      <,    '  '  '  ' 

HIDROLOGIC  INSTBUIVIENTATION  ■""  -      '    "    '■-"'V-:--  '    -  ■:■ 

During  the  summer  of  1963  recgrding^  and  non-reeording  rainga^es  were 
installed  at  selected  ranches  on  each  watershed^  and  water  level  recorders 
were  placed  in  the  stremn  near  each  USGS  crest-stage  gagec    All  the  raingages 
were  built  to  UcSo  Weather  Bureau  specifications j,  and  were  installed  in 
accordance  with  their  suggestions o    Shelters  and  stilling  wells  for  the  • 
water  level  recorders  were  made  from  designs  adapted  from  UoSoGoSo  drawings „ 
Details  of  design  md  installation  were  included  in  Interim  Report  #2o 


(Williams  J,  19651b  )o    A  view  of  one  of  the  water  level  recorder  installations 
is  shown  in  Figure  26  o 

Weather  station  instrument  packages  were  developed,  by  the  Electronics 
Research  Laboratory  at  Montana  State  University^  and  eight    of  these  iinits  were 
installed  at  ranches  on  the  watersheds  during  the  e\iiimier  of  196^*0  Each 
package  consists  of  a  group  of  sensors  to  measure  wind  speed,  wind  direction, 
air  temperatures  soil  temperature  and  soil  moisture;  semi-automatic  recorder 
to  record  a  signal  from  each  sensor  once  every  30  minutes;  and  cable  strung 
underground  to  connect  the  sensors  to  the  recordero 

.  At  each  ranch  the  sensors  were  located  in  a  fenced  enclosure  10  ft  X  10 
ft  square  5  situated  between  200  and  500  feet  from  the  ranch  house  <> 

Wind  speed  and  direction  were  measured  by  a  cup  anemometer  and  single-tail 
direction  vane  built  by  the  Belfort  Instrument  Company  c    The  wind  instrments 
were  mounted  atop  a  pole  16  feet  above  the  grounds    Air  and  soil  temperatures 
were  measured  by  thermistors  made  by  the  Yellow  Springs  Instrument  Company o  ' 
The  air  thermistors  were  mounted  between  two  aluainim  discs  to  permit  free 
air  movement  and  to  obtain  the  proper  radiation  environment ^  while  the  soil 
thermistors  were  encapsulated  in  epo3cy  for  protection »    Soil  moisture  was 
measured  by  ^f^sxm  soil  blocks  until  196Tj  and  after  that  time  by  plastic 
Boijyoucos  blocks o    Either  style  of  block  contains  two  electrodes  to  which 
alternating  voltage  is  applied ^  the  resulting  current  being  a  function 
©f  the  moisture  in  the  surrounding  soilo 

Considerable  electronic  circuitry  was  required  to  convert  the 
outputs  from  the  various  sensors  into  a  form  suitable  for  recording »  The 
electronic  equipment  and  single-channel  Rustrak  recorder  were  housed  in  a 

„  1+8  - 


FIGURE  26  —  Water  Level  Recorder  Installation. 
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metal  cabinet  about  l8  inches  by  10  inches  by  10  inches  o    The  cabinet  was 
placed  in  the  rancher" s  houses  usually  in  a  basement  or  storeroomc  The 
system  was  designed  to  provide  one  recording  from  each  of  ten  sensors  at 
30-minute  intervals c    An  interval  timer  was  used  to  start  and  stop  the 
station  automatically  while  a  gold-plated  telephone  type  stepper  switch  was 
provided  to  automatically  switch  (commutate)  the  output  from  each  sensor  in 
turn  into  the  recordero 

Output  frc»  the  Rustrak  recorder  is  a  2o 5-inch-wide  pressure  sensitive 
paper  tapgo    Since  the  recorder  was  only  in  operation  for  a  short  time  every 
30  minutes g  a  roll  of  paper  tape  lasted  for  one  month o    The  outputs  from  ten 
sensors  appeared  as  a  series  of  short  marks  on  the  tape,  with  about  one 
inch  of  tape  "being  used  every  30  minutes  o  Identification  of  the  outputs  on 
the  tape  was  possible  because  the  wind  speed  which  was  always;  the  first 
channel  of  the  output  ^  produced  a  slanting  mark  across  the  tape  2,  while  the 
other  instruments  all  produced  longitudinal  marks  and  always  recorded  in  the 
same  seq^ueneeo    A  recording  timer  in  the  circuitry  caused  the  recording  time 
to  be  doubled  once  every  2^  hours j  so  the  midnight  recording  appeared  as  a 
series  of  elongated  markes c    A  sample  of  recorder  output ^  and  the  calibration 
codes  are  shown  in  Figure  27  <> 

Details  of  the  design  and  installation  of  the  weather  station  packages 
were  given  in  Interim  Report  #3  (Williams  and  Edwards ^  1965)° 

Figures  28  and  29  show  one  of  the  weather  station  sensor  enclosures g 

and  one  of  the  recorder  units « 
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Several  operational  problems  arose  with  the  reeording  equipment j  and 
in  1967  the  Electronics  Research  Laboratory  made  extensive  modifications  to 
improve  reliability  and  accuracy o    Because  additional  components  were  added^ 
larger  cabinets ^  2^^  inches  by  12  inches  by  I8  inches,  were  requiredo 

Figure  30  shows  one  of  the  rebuilt  recorder  units o 

As  originally  installed  only  four  of  the  weather  stations  were  equipped 
to  measure  wind  speed  and  direction  =>    In  196T  wind  equipment  was  added  to 
the  four  remaining  stations o 

Data  Reeording 

Arrangements  were  made  with  each  of  the  ranchers  where  instruments  were 
installedj  for  them  to  change  charts  on  the  instrmentSj  wind  clocks  ^  etCo 
Non«=re cording  raingages  required  daily  inspection  and  measurement  of 
contents;  reeording  raingages  required  weekly  inspection ^  clock  winding  and 
chart  removals  and  occasional  emptying  of  the  bucket;  weather  stations 
required  monthly  chart  removal,  and  periodic  inspection  to  see  that  the 
instrument  was  operating  satisfactorilyo    Ranchers  were  paid  a  nominal 
fee  annually  for  their  services o 

A  cooperative  agreement  was  made  with  the  UoSo  Geological  Survey 
whereby  the  USGS  engineers  regularly  inspected  the  water  level  recorders  5, 
wound  clocks  9  changed  charts  and  performed  maintenance  service  as  requiredo 
The  USGS  visits  coincided  with  their  inspection  of  the  crest-stage  gages 
so  that  extra  trips  were  not  necessary 0    The  USGS  also  made  current 
meter  measurements  when  possible  to  improve  their  established  stage -discharge 
relationships o    Cheek  levels  were  run  by  the  USGS  to  each  recorder  instal- 
lation annually  o 
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FIGURE  30  —  Rebuilt  Weather  Station  Recorder. 
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Table  II  shows  the  extent  of  installation  on  each  of  the  watersheds o 


TABLE  II 


Area  in 

Reeordinc 

Non— recording 

Weather 

Water  Level 

sq^o  mi 0 

Rain  Ga^es 

Rain  Gages 

Stations* 

Re e orders 

Baeon  Creek 

I80O 

2 

3 

2 

1 

Duck  Creek 

53o8 

3 

k 

3 

2 

Hmp  Creek 

7o6 

1 

1 

1 

1 

Lone  Man  Coialee 

Ik. 6 

2 

■  2 

2 

1 

*Each  weather  station  equipped  to  measure  soil  temperatijsre  and  soil 
moisture  at  three  soil  depths ;  air  temperature  at  two  elevations  9  wind 
speed  and  wind  direction o  (Wind  measurements  were  omitted  from  k  of  the 
stations  when  originally  installed;  wind  was  added  to  these  four  stations 
in  1967) 


A  "brief  summary  of  the  hydrologic  data  gathered  is  presented  in 
Appendi2c  Eo  >-"-;      .  ■ 

SNOW  STUDIES 

Two  techniques  were  employed  for  evaluating  snow  melt  potential 
on  the  watersheds o    These  were  periodic  snow  surveys,  and  aerial  photos o 

The  State  Snow  Survey  Supervisors,  UoSo  Soil  Conservation  Service  aided 
in  the  establishment  of  three  snow  survey  stations  on  each  of  the  four 
watersheds  o    SCS  personnel  from  the  local  work  units  then  made  surveys 
within  a  few  days  after  each  snowstorm  of  six  inches  or  moreo    At  each 
of  the  three  locations  the  surveyor  collected  10  samples  with  a  two-inch 
core  smpling  tube  5  recording  snow  depths  and  core  weights  o    Th®  core  weights 
were  then  converted  into  inches  of  water o 
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Smma,Tj  snow  iTurvey  data  are  shown  in  Appendix  E    o        '■  ' 

Aerial  photographs  were  taken  by  the  Montana  Highway  Department  photo- 
graphy  seetion  from  the  department's  Aero-Comandero    It  had  been  hoped 
that  flights  eould  be  made  after  eaeh  major  storm  to  evaluate  the  effeet  of 
driftingo    Weather  conditions  and  scheduling  problems  prevented  the  obtaining 
of  photographs  as  frequently  as  desired^  but  one  set  vae  obtained  for  eaeh 
watershed  in  February  or  Mareh  of  196h^  I965,  1966 ^  1967  and  19680    These  ■ 
photos  were  used  qualitatively  to  estimate  percentage  of  watershed  that  was 
snow  coveredo    Althou^  the  quality  of  the  photos  would  permit  stereoscopie 
viewingg  no  estimates  of  snow  depths  were  madec  '  ■■  . 

ASSE3SME1T  OF  EFFECTI¥ENESS  OF  HYDROLOGIC  DATA  COLLECTION  'r^  /::.  ■.^^i'V/ ^-V:. 

Precipitation  and  streamflow  data  were  collected  on  the  four  water- 
sheds for  si3c  years g  from  September  19 6 3  to  September  1969°  Supplemental 
data  were  collected  on  three  of  the  watersheds  .from  April  to  June  1970 o'  ■ 
Additional  meteorological  information  (wind  speed  aad  directions  air  and  soil 
temperature 9  soil  moisture)  was  obtained  between  July  196k  and  September  ' 
19690    Seme  comments  regarding  effectiveness  of  the  data  collection  program 
follow  o  .      •       ■.  . 

Precipitation  Data  '  ■,  ■  ■ 

There  were  relatively  few  gaps  in  data  at  most  stations  during  the 

six  years o    One  rancher- observer  at  Hump  Creek g  where  a  recording  raingage 

was  establishedg  was  frequently  absent  from  the  ranch  for  extended  periods  9 

and  the  8°=day  clock  on  his  raingage  was  often  allowed    to  run  downo  This 

station  had  more  missing  records  than  any  other o  ,■ . 
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Reliability  of  the  precipitation  data  is  felt  t©  be  generally  y@tj  good 
when  the  precipitation  fell  as  rain 9  and  fair  to  poor  when  it  fell  as  snowo 
Snow  on  these  watersheds  is  almost  always  aecompanied  "by  wind  of  varying 
intensities  J  md  the  precipitation  caiight  by  the  raingages  under  such  eon<" 
ditions  is  often  ae  low  as  50  percent  of  the  true  precipitation o 

The  efficiency  of  the  raingages  during  winter  months  could  perhaps  have 
been  improved  by  the  addition  of  windshields o    The  decision  t©  omit  wind-' 
shields  was  based  on  economics  ^  and  the  resTult  of  conversations  with  several 
researchers  at  MsU  who  generally  reported  that  the  use  of  windshields  had 
apparently  feeea  only  partially  effective  in  their  studies o 

Streamflow  Data         .  ■ ,       .  . 

Streamflow  records  obtained  for  the  four  watersheds  are  believed  to  fee  very 
goodo    The  only  major  missing  data  is  from  the  Duck  Creek  station  for 

The  Duck  Creek  stilling  well  and  recorder  were  originally  attached 
to  a  county  bridge  over  Duck  Greeks  but  the  bridge  was  so  light  that  traffic 
caused  serious  vibrations  and  made  recorder  operation  unfeasableo    In  196^ 
the  recorder  was  relocated  away  from  the  bridge  and  no  further  difficulties 
were  encounteredo 

During  winter  months  the  float  in  each  of  the  stilling  wells  was 
suspended  out  of  the  water  to  prevent  freezing ^  in  such  a  way  as  to  float  ^ 
free  when  the  water  level  rose;  it  is  believed  only  very  minor  streamflow 
events  were  missedo 

There  were  long  periods  at  each  of  the  stations  with  no  floWg  or  almost 
no  floWo    Gaging  operations  might  have  "been  more  efficient  had  modifications 
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been  made  to  the  reeorders  to  prevent  them  from  manning  Tuntll  the  water 
level  roseo    An  attachment  for  this  purpose  is  available  g  and  gaiises  a  .,,    :  , 
reeorder  to  sit  idle  until  streraiflow  oeeurSg  at  whieh  time  it  eoamenees  .... 
operation  and  rms  naatil  an  inspector  stops  it  after  the  stresmflow  eeasesc 

Weather  Station  Packaiges 

The  weather  station  paekages  were  developed  at  Montana  State  University 
because  comparable  equipment  was  not  available  commercially  at  prices  the 
Project  could  af fordo    The  eqiuipaient  as  developed  provided  a  ma^imiM         ■  ■■: 
amoTumt  of  meteorological  data  for  the  price;  however,  eonspromises  were  ■  ■-.  .- 
necessary 9  and  severa\l  problems  arose  as  a  result  o    Althou^  operation  was 
improved  when  the  electronic  packages  were  rebuilt  in  19^7^  bomb  of  the 
problems  never  were  eliminated  entireljo  . 

Changing  charts  on  the  RustrA  recorder  is  an  intricate  process  which 
mTiist  be  done  very  carefully  if  the  recorder  is  to  operate  properly o  The 
chart  drive  operates  through  a  friction  clutch  j  and  unless  the  char's  • 
is  installed  perfectly »  it  will  not  advance o    The  rancher-observers  on  the 
project  had  considerable  difficulty  with  the  charts  o  • .  ; 

Wien  the  umits  were  rebuilt  in  196T  the  recorders  were  installed  so 
they  could  be  removed  by  'onpluggingo    Two  spare  recorders  were  purchased  ■ 
so  a  rimeher  could  unplug  oneg  mail  it  to  the  project  supervisor g  and  install 
smother  recorder  on  which  the  diarts  was  already  installedo    For  a  time  the 
recorders  were  mailed  back  and  forth  regularly ^  but  the  operation  did  not 
prove  to  be  much  more  reliable  even  with  "experts"  installing  the  charts o 


0<seasi©nal  power  fadlures  at  the  ranches  eaused.  the  etationg  to  fee  out 
of  operation  for  varying  periods  o    This  eaiaeed  the  "midnight  elongation"  of 
the  output  to  be  generated  at  the  wrong  time  ^  mitil  the  ©"b server  re-'Set  a 
timer  on  the  mito    This  problem  was  allevJ.ated  in  1967  when  the  timer  was 
converted  to  flashlight  battery  powero 

With  the  exeeption  of  soil  moistwe  elements  ^  very  few  problems  were 
eneountered  with  the  sensors  situated  in  the  instrument  enelosures  o  Both 
t^fpes  ©f  soil  moisture  block  used  (gypsum  and  Bouyoucos)  operate  by  measuring 
current  between  two  electrodes  embedded  in  the  element  in  the  presenee  of 
alternating  voltage o    The  current  varies  with  moisture  tension  in  the 
soil  3,  and  the  tension  is  in  turn  related  to  the  moisture  content  o    The  blocks 
were  calibrated  to  establish  a  current -tens ion  relationship  for  eaehg 
and  soil  samples  from  each  of  the  installation  sites  were  tested  to 
determine  tension-moisture  relationships  o    The  ^^sum  blocks  ^  which  were 
originally  installedg,  were  thought  be  be  superior  for  use  in  some  Montana 
soils  which  are  high  in  salt  content;'  some  of  these  blocks  deteriorated 
rapidly  (probably  owing  to  freeze'-thaw  conditions)  so  they  were  all 
replaced  in  196J  by  the  more  durable  Boti^oueos  blocks  0 

Calibration  of  the  soil  moisture  elements  was  never  entirely  sat- 
isfactory; soil  moisture  data^  therefore  9  should  be  used  qualitatively  only 9 
to  establish  periods  of  relatively  wet  or  relatively  dry  soilj  and  not  to 
evaluate  precisely  the  soil  moisture  conditions  c 

The  output  records  from  the  RustrA  recorders  were  converted  into 
digital  form  by  manual=- visual  means  o    A  minimm  of  8  man»»hours  of  student 
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labor  was  regmred  for  the  reduetlon  of  data  from  one  chart  o  Inea^erieneed 
students  took  up  to  twice  that  momit  of  timeo 

SUMMARY 

The  procedures  followed  in  the  hydrologic  studies  of  four  small 
■watersheds  hare  been  deseribedo    Data  collection  and  reduction  probleais 
have  been  discussedo    Siaamaries  of  watershed  characteristics  t  aad  hydrologi 
data  gather  are  presented  in  Appendices  Bj       D  and  Eo 
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CHAPTER  lY 


THE  19 61*  FLOOD  AT  LONE  MM  COULEE 


On  Jime  196k  a  rain  stom  of  impreeedented  magnitude  and  areal 

extent  struek.  a  wide  area  of  north-eentral  and  northwestern  Montaaao  The 
stom,  which  was  centered  over  the  continental  divide  g  resulted  in  floods 
nruieh  larger  than  my  which  had  ever  been  recorded  on  umj  watersheds  o  Lone 
Man  Coulee  J,  in  Pondera  County  is  located  about  50  miles  east  of  the  storm 
centers  but  it  received  over  5  inches  of  rain  during  the  30-hour  stonaao  The 
resulting  hydrograph  had  a  pe^  flow  of  17^0  efs  md  total  runoff  of  I067 
inches o    The  watershed  had  been  selected  in  1963  as  one  of  the  Project 
watersheds  to  be  studied  comprehensively  in  this  investigation j  and  during 
the  196^*  stom  four  radngages  (two  recording  and  two  non-re cording)  and 
one  continuous  recording  water  level  recorder  were  in  operation o         ;.  • 

The  Project  supervisor  visited  the  watershed  on  June  89  196k ^  and  took 
a  number  of  pictures  at  the  mouth  of  the  stream  shortly  before  the  failure  of  a 
six-foot  culvert  under  a  county  roado    See  Figures  31  to  3^o 

Project  personnel  returned  to  the  watershed  a  week  after  the  floodg 
flagged  high-water  marks  g  and  conducted  a  slope-area  study  to  determine 
discharge  peako 

The  hydrology  of  Lone  Man  Coulee  was  studied  in  detaa,l  using  data  from 
the  June  196^  flood  and  other  runoff  events  which  ©ccarred  in  196k  and  ,1965 0  TV© 
reports  of  fin^dings  were  prepared  (Williams  ^  196^  and  Williams  ^  Bobinsoa 
md  Hanson  5,  1966)0 


Lone  lyian  Coulee  ¥S  Recorder  and  Culvert  Fill(  Note 
Man  Crossing  Fill  )  6/  8/64 


FIGUEE  31» 
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Lone  Man  Coulee  -  Across  Channel  from  Left  Bank 

6/l8/6ii. 


Downstream  ©veraH  "View  of  Lone  Man  Reach 

6/18/64 

FIGURE  32. 


Lone  Man  Coulee  -  Looking  Upstream  from  Blirff  above 
Culvert  (  Mller  Coulee  from  Right  )  6/  Q/6k 


Upstream  at  Lone  Man  Coulee  (  Confluence  with  MiHer 
Hidden  by  Rocks  ) 


FIGURE  33 
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Lone  Man  Coulee  -  Flooding  Culvert  Upstream  from 
¥S  Recorder  and  Culvert  Fill  6/  8/6^ 


Lone  Man  Coulee  -  Looking  Downstream  at  WS  Recorder 
gnd  Culvert  Fill  6/l8/6i^ 


FIGURE  3^ 
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CHAPTER  ¥  -     -  ■  • 

MULTIfAEIATE  STATISTICAL  STUDIES 

Prineipal  (somponent  analysis  and  varlmax  rotation  of  the  principal 
factors 5,  ■yising  watershed  >gharaeteristic®  and  bydrologie  data  determined  for 
the  four  projeet  watersheds ^  were  performed  to  provide  information  aboMt  the 
relative  import anee  of  29  independent  Tariables  to  the  peak  diiseharge  rates 
and  nmoff  volimes  produced  hj  these  variables  o    A  aeries  of  eonapnter  pro-gr^is 
was  written  to  perform  the  TarioMB  operations  o    Data  from  50  rmoff  events 
were  'used  9  twenty  of  whieh  were  snowmelt ,  and  30  were  rain-ea^agedo         ■•  ■■■ 

The  procuduras  are  described  "briefly  in  this  ehapterg,  and  the 
results  are  diseussedo    The  proeediores  and  the  reasons  for  their  seleetioUg 
together  with  theoretical  'baekgroiand  material  were  presented  in  detail  in  ■■-  ■ 
Interim  Report  #6  (Lewis  and  Williams  ^  1968)0    Since  piib ligation  of  that 
report  further  aaalysis  has  "been  made  laeing  the  same  data^  and  a  e«pMter 
error  whieh  eauses  slight  changes  in  the  resToltSg,  has  "been  eorreetedo    ;  , 

REVIEW  OF  PROCEBUBES  v  —  '  ' 

Five  watersheds  were  considered  (East  Duck  Creek  aad  Budk  Creek  were 
treated  as  separate  watersheds  because  they  were  each  eqiaipped  with  a  . water 
level  recorder)  o    Twelve  watershed  parameters  each  of  which  were  assuaed  t© 
remain  constant  with  timeg  were  determined  for  each  Casino    Field  meaenre'- 
mentSg  as  well  as  map  and  aerial  photograph  analysis  were  eia^loyedo    In  some 
cases  weighted  mem  values  were  lasedo  ■       '  .■.  - 
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SeTenteen  st©m  Tariables  were  deteiml.ned  for  eaeh  of  fifty  rmoff 
events  whieh  oeearred  oa  the  five  watersheds  duaring  the  years  1964=^67° 

Deseriptions  ©f  the  29  parameters  aad  tabialated  vsilmes  for  eaeh  are 
shown  in  Appendix    F  o       ...  .  ■  -  . 

The  amomt  aad  type  of  data  aTailatole  sm^ested  that  a  prineipal  ©om- 
ponent  extraetioa^  followed  "by  a  Tarimax  rotation  for  interpretstioa  of  the 
variables g  aad  a  mltiple  regression  on  either  the  prineipal  eomponents  or 
the  rotated  faetors  for  the  predietion  equation  would  represent  the  Ibest 
statistical  system  of  analysis  o     (See  Lewis  and  WillimSa  19685,  for  eoMents 
of  various  infestigators  proposing  these  methods).,  - 

The  variolas  steps  repaired  in  the  stuady  are  shown  in  bloek  diagrua  form 
in  Figiore  35  ^  aad  are  described  below o    Eaeh  step  is  a  eompinter  prograsa 
solved  sequentially g  each  Btiliiing  the  resets  of  preceding  steps  o 

Correlation  is  a  program  to  eoi^ute  means  ^  standard  dematioas^  wr-- 
iances  and  eo^variances §  Bvms  ©f  squaares 5  and  correlation  coefficients  of 
any  n«ber  of  observations  for       to  32  different  variables  or  parameters  o 
(See  Lewis  &  Willi ms,,  1968^  Appendix  S)o 

Principal  Coasgonent^^  Analysis  is  a  program  which  perfoms  a  principal 
component  aaaaljsis  on  a  correlation  matrix  (mp  to  60  1  60)  mid  prodiaces 
Top  to  60  nomaliied  eigenvectors  and  the  respective  eigen"vali!aes  of  each 
factoro    The  principal  factors  are  also  computedo    The  variance  (eigenvalue) 
of  each  vectors  percent  of  variance  e^lained  by  each  factor ^  and  aceraafiilated 
percent  are  also  compnitedo     (See  Lewis  &  Williams  5,  19685,  Appendix  F)  o . 


Raw  Data 


Correlation:        ,  , 
Means , 

Standard  Deviations, 
Logarithms 


Principal  Components 
Analysis 


Regression  of 
Principal  Components 


Varimax  Rotation 


Output 


Standardization 


Factor  Value 
Computation 


One-Pass 

Multiple  Regression 


Regression 
Coefficients  for 
Variables 


FIGURE  35  --  Block  Diagram  of  Multivariate 
Procedures 
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RegrgsglQn^^  of  Pginelpal  CQwponentB_  is  a  progrrai  i&ieh  eqmp^tes 
prineipal-eaapoaent  regression  eoefficients  for  the  eigenvectors o 

Yarimay  Rotatiosi  is  a  program  to  perform  an  orthogonal  normal  varim^ 
rotation  ©f  as  msnj  m  75  faetorg  with  loadings  for  as  maaj  m  100  varia'bles  o 
Variances  J  pereents  and  aeemulated  pereents  of  the  final  rotated  faetors 
are  aleo  eompiatado     (See  Lewis  and  Williams  g  1968 5  Appendix  G) 

Standardisation^  is  a  prograM  whieh  eomputes  the  stMidardiied  form  of  31 
variables  from  the  observations  5,  variable  means  and  variable  standard 
deviations  o    The  new  values  are  dimensionless ,  hm®  a  Man  ©f  lerog  and 
a  standard  deviation  of  imityo  .    •  ■  ■ 

Faetor  Ya-l^e  Go^gtation  is  a  program  whieh  computes  the  values  of  1^ 
to  10  factors  lasing  data  from  up  to  50  obserrations  of  as  maaj  m  31 
variables o    The  operator  seleets  whieh  of  the  31  variables  are  important  mid 
are  to  be  med  in  the  f aetor  equations  g  by  st'adying  the  ©latpiit  from  the 
varimax  rotation  programo  :■  ,  

One-^Pass  Multiple  Regression^  is  a  program  which  computes  means  ^ 
variances  a  stmdard  deviations  and  all  pairs  of  correlations  for  as  many  as 
30  variables o    Any  one  of  the  given  set  can  be  selected  as  the  dependent 
variable  s  with  any  sub^set  of  the  remaining  variables  used  as  the  independent 
variables  o 
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Bg£rgggioia  CQ©ffloien;^s  for  Variables ^  is  a  program  which  eonverts  the 
regression  eoeffieienti  for  faetors  into  eoeffieients  for  the  variables 
&©utm.ned  in  the  faetors  <= 

RESULTS  •• 

The  principal  coisponent  analysis  resulted  in  the  generation  of  exponents 
B^j,  B^g  By  etoo  in  the  eqiaations 

and  the  constant  B  o    In  this  equation  I  is  the  dependent  variable  (peak 

o 

discharge  rate  or  total  runoff  voliMe)  and  X^,  X^j,  X^,  etCo  are  values  of  the 
29  watershed  parameters  and  storm  variables  o    Table  III  (a)  preeents  the 
applicable  ea^onents  and  constant  for  the  analysis  using  all  50  events g  and 
al^o  for  a  later  analysis  ijsing  the  30  events  that  were  rain-eainsedo  (With 
the  30  events  snow  water  equivalent  was  omitted s  leaving  28  variables )o 
In  e^^ining  T^le  III  it  should  be  noted  that  negative  signs  iicply  inverse 
correlation;  io@os  peak  discharge  rate  varies  inversely  with  FREQ,  (stream 
fret^ency)  and  directly  with  ELE¥  (watershed  elevation )o    The  principal 
coiiponent  analysis  does  not  provide  information  as  to  the  relative  importance 
©f  each  independent  variable  o       ^;     - -■ 

The  purpose  of  the  varimax  rotation  was  to  "rotate  each  of  the  principal 
cQSiponents"  s©  that  each  component  woiO-d  have  high  coefficients  on  certain'; 
variables  and  low  coefficients  on  the  other  variables  5  thereby  allowing 
interpretations  for  the  important  variables o    Using  data  from  all  50 
events 3,  varimax  rotation  classified  the  29  independent  variables  into  ten  , 
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groups  ©r  "factors"  and  indicated  which  variables  were  relatively  more  - 
importaioit  to  eaeh  faetor^    Also,  each  factor  is  relatively  less  important 
thm  the  preceding  factors  o    The  same  procediire  was  followed  for  the  30 
raiB-eaiased  events s,  classifying  the  28  independent  variables  into  11  factors o 
The  varimax  rotation  results  are  shown  in  Table  IV  «    The  variables  deemed 
to  be  most  important  are  shown  in  the  shaded  areas  of  the  table »  .;: 

Following  varimax  rotation  and  perusal  of  Table  IV  g  any  desired 
combination  of  variables  can  be  used  in  factor  value  computation  and 
regression  analysis  <,    This  was  done  with  l8  selected  variables  and  all  50 
events  9  md  again  with  the  30  rain-caused  events  o    As  with  principal 
eeaponent  analysis  the  result  is  generation  of  esqponents  B^^  B^g  B^s,  etco 
and  the  constant  B^  in  equation  8.    Table  III '  ("b')' 'Shows  the  meults  of 
this  analysis  o    As  before  ^  the  sign  of  each  exponent  is  more  significant 
thas.  its  magnitude o  '       ■  ."■    ■      ■  -  '     ''  '  ' 

After  publication  of  Interim  Report  #6  (Lewis  and  Williams g  I966)  a 
computer  progrOTiming  error  in  the  factor  value  computation  program  was 
detected^  and  this  is  responsible  for  slight  errors  in  Table  XI  (page  7^) 

of  that  report,  c  -   ■  ■ 

DISCUS  SI  01  OF  RESULTS  .       .       ,.,  ,  .„  ,  ;.,:.:■.„ 

A  search  of  literature  failed  to  show  that  the  techniq.ues  used 
in  this  phase  of  the  investigation  had  been  used  before  for  hydrologic 
studies  as  large  as  this  oneo    These  techniques  are  much  more  useful  than 
ordinary  regression  methods  for  two  reasons o    Firsts  it  is  presupposed  in  a 
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TABLE  IV:     Successive  Importance  of  Variables  to  Factors 


a)    Fifty  Events 


Factor 


Decreasing  Order  of  Importance 


b)    Thirty  Rain-Caused  Events 


Fa.  c  tor 


Decreasing  Order  of  Importance 


Note:     The  variables  deemed  to  be  most  important  are  shown  in  the 
shaded  areas  of  the  table. 

Parameters  are  identified  in  Appendix  F. 
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regression  studj  tliat  all  the  variables  used  are  mutually  independent »  In 
a  hy'drologic  study  this  is  not  a  valid  assumption  because  it  cannot  "be  stated 
for  inatanee  that  precipitation  intensity  oeciars  eompletely  independently  of 
wind  veloeityo    To  the  extent  that  the  variables  are  inter- relatedj  ordinary 
regression  methods  will  give  erroneous  results »    The  techniques  selected  for 
this  investigation  have  "been  shown  mathematically  to  eliminate  the  effects 
of  inter=»eorrelation  among  the  v&riables  o    Second,  ordinary  techniques  - 
cajanot  indicate  the  relative  importance  of  the  variables  usedo    The  varimax 
rotation  with  factor  analysis  employed  in  this  study  provides  a  means  of 
determining  which  variables  are  most  important o  •..  .  ■ 


Multivariate  analyses  were  used  in  this  investigation  to  see  whether  •  -. 
they  offered  possibilities  for  delineating  the  important  parameters  affecting 
flood  peak  magnitudes,  as  well  as  to  determine  the  form  of  the  regression 
equation  for  flood  pe^  magnitude o    Results  of  the  study  must  be  considered 
preliminarja  but  the  objectives  were  mete    The  study  showed  that  it  is  indeed 
possible  in  this  manner  to  delineate  important  parameters  ^  and  it  is  possible 
to  develop  a  regression  equation  for  flood  peak  magnitudeo     ■  ■     ■     ■   -.  ■  .  ■ 

There  are  several  ways  in  which  the  study  could  be  improved  i^on^  if  it 
were  to  be  done  againc    These  includes    ,■  - 

lo    Redefine  some  of  the  variables o    Precipitation  intensity  should 
be  defined  m  the  average  intensity  during  the  period  of  excess  rainfall  of 
the  storm  instead  of  the  average  intensity  computed  for  the  beginning 
©f  rainfall  to  end  of  nmoffo    Stom  duration  j  defined  as  average  time  to 
center  of  area  of  all  rainfall  hydrographs  for  the  watershed  does  not  appear 
t©  be  meaningful  in  terms  ©f  peak  discharge  magnitude  j  althoiagh  it  may  •: 
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"be  related  to  total  nmoff  volxaneo    Total  precipitation  likewise  is  probably 
not  meaulngful  in  establishing  peak  discharge  rates o 

■  2o    Provide  a  more  dense  network  of  instrimentationo    Sane  parameters 
were  based  t^on  single  or  a  fev  measurements  even  thou^  they  were  quite 
variable  aeross  the  watershedo 

■. .  3o    Include  geographical  factors  o    With  the  exception  of  watershed 
elevation  no  parmeters  describing  geographical  location  within  the  state 
were  inclmdedo    '         -  ■  - 

ko    Try  more  combinations  of  parameters  in  the  factor  analysis  c  There 
mgQT  be  a  combination  of  variables  that  is  better  than  the  l8  used  in  the 
preparation  of  Table  III  (b)  (or  the  ih  and  6  reported  by  Lewis  and  Williams  9 
1968)0  -    -  -    -  ■  . 

Consideration  of  the  es^onents  resulting  from  principal  component  as 
compared  to  factor  regression  on  Table  III  shows  that  the  maximum  discharge 
^tiadj  of  30  events  s  the  two  methods  result  in  opposite  sign  of  e^onent  on 
five  variables  3,  L  (stream  length)  ^  USE  (land  use  ratio)  5,  API  (antecedent 
precipitation  index) 5  WEEK  (week  of  the  year)  and  AIRT  (air  temperature ) o 
Intuitively g  the  sign  for  API  seems  correct  on  principal  component  regression 
and  wrong  ©n  factor  regressiono     (The  positive  sign  indicates  increasing 
peak  discharge  rate  with  increasing  magnitude  of  APl)o    It  is  not  immedi- 
ately  evident  that  either  sign  is  better  for  the  other  four  exponents c  Two 
exponents  appear  intuitively  wrong  for  both  methods  GNDS  (overland  ground 
slope)  Bnd  S  (  main  channel  slope )o    Discharge  rate  would  logically  vary 
directly  with  both  these  terais  o 


Examination  of  the  listing  of  relative  in^ortance  of  variables  to 
factors  (Table  IV  )  shows  that  the  study  has  identified  the  obviously  import- 
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ant  parameters  for  the  analysis  of  30  eTentSo    When  50  events  are  considered^ 
the  effect  of  lumping  snowmelt  and  rain-caused  events  together  has 
masked  the  importance  of  either  SOLM  (soil  moisture)  or  API  (antecedent 
precipitation  index)  o    Certainly  with  rain-caused  events  one  of  these  vari=- 
ables  should  be  ineludedo       •-■■^        ■  '        •  •  •■  ^- 

Interim  Report  #6  (lewis  and  Williams,  I968)  describing  the  multivariate 
statistical  studies  was  reviewed  by  Do  C  Woo,  for  the  Structural  and  Applied 
Mechanics  Divisions,  Federal  Highway  Administration o    Woo  had  comments  on  two 
araass    input  data^  and  hydrological  interpretationo ' 

Regarding  input  dataj  Woo  noted  that  no  assessment  of  the  acciiracy  of  the 
recorded  data  was  attempted.    Future  studies  of  this  nature  g-hould  include 
validity  checks  on  the  data  collectedo    This  willg  however^  require  a  greatly 
expanded  ^ata  collection  network j  and  hence  a  more  costly  study o 

Regarding  hydrological  interpretation.  Woo  suggested  that  runoff  events 
should  be  groined  by  magnitude g  and  events  of  comparable  magnitude  considered 
separately o    He  further  pointed  out  that  snowmelt  events  should  be  treated 
apart  from  rain<=<sauBed  events  o    Future  studies  should  consider  classification 
of  events  by  magnitude  as  Woo  suggests  o    As  reported  in  this  chapters,  rain- 
caused  events  have  been  studied  separately  since  publication  of  Interim  Report 
#60 

CONCLUSIONS  ;.. 

The  results  of  the  analysis  described  here  point  to  a  recommendation 
that  in  the  future  hydrologic  study  of  small  watersheds,  where  peak  discharge 
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detemintation  is  the  goal,  attention  should  be  focused  to  obtaining  as  much 
reliable  information  as  possible  on  precipitation  intensities  and  amounts  5 
and  air  temperatures;  with  lessened  emphasis  on  wind  speeds  and  directions 9 
and  soil  moisture  measurements »    Soil  temperature  measurements  are  important 

when  snowmelt  is  considered,  but  unnecessary  when  rain-=caused  events  only 

are  considered,. 

,■  A  further  reconmendation  is  that  principal  component  analysis  be  tried 
with  a  lesser  number  of  parameters.,    The  results  of  this  study  verify  the 
contention  of  Amoroeho  and  Hart  (196^1)  that  the  use  of  fewer  parameters  is 
generally  more  successful  than  the  use  of  many  parameters c  Reduction 
of  the  number  of  parameters  would  rule  out  the  use  of  factor  analysis j 
howeverg  because  20  is  considered  the  smallest  n^lmber  which  will  produce 
satisfactory  results  in  such  a  study o        <  v-  ■  ■  I :;  ■  ... 

This  chapter  has  described  the  procedures  used  in  multivariate 
statistical  analysis  of  data  from  five  watersheds  and  fifty  runoff  events ^ 
and  has  discussed  the  results  obtainedo 

"■' ...       •  ■ 
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•i^^p.j  -^:     •  :•       ;:>.'    CHAPTER  VI 

'       STUDY  OF  SOIL  CONSERVATI.QII  SERVICE  METHOD 

The  UcSo  Soil  Conservation  Service  has  developed  several  methods  . 
for  determining  peak  flow  rates  from  small  watersheds.    Some  of  these 
methods  are  discussed  in  Chapter  VII  of  this  report;  but  the  principal 
method  currently  in  use  by  the  SCS  is  reviewed  in  this  chapter,  and  the 
results  of  a  test  of  the  methodj  using  data  from  an  actual  storm  which 
occurred  on  Duck  Creek  watershed,  are  reportedo    More  details  are  supplied 
in  Interim  Report  #5  (Ferris  and  Williams,  I968),  .  .  v 

SCS,  METHOD        ;;;;  .--■^^  ,  ■^■j- 

In  their  watershed  planning  studies.  Soil  Conservation  Service  . 
engineers  are  first  interested  in  total  runoff  volume,  because  total 
volme  determinations  are  critical  in  establishing  the  feasibility  of 
reservoir  construction  and  other  reclamation  activities o    Discharge  rate  is 
of  importance  only  insofar  as  it  affects  spillway  capacities,  etCo    In  all 
their  hydrologic  methods,  therefore,  determination  of  total  runoff  volume 
is  accomplished  first ,  using  methods  developed  by  SCS  engineers «  Various 
methods  are  then  employed  for  determination  of  peak  discharge  rate;  these 
methods  usually  involve  standard  flood  routing  techniques ,  unit  hydrograph 
construction g  or  others  o  ^'i-j „:>;;;     "-.i'..'--      "  .v;  -  '''^  ■:  .  : 

None  of  the  SCS  methods  assign  a  return  frequency  to  the  runoff  volume 
or  peak  flow  rateo    Rather,  they  work  from  a  real  or  hypothetical  rainstorm 


-  81  - 


and  although  the  rainstorm  may  "be  assigned  a  return  frequency,  the  resulting 

runoff  event  is  not., 


Runoff  volumes 

The  accumulated  volume  of  runoff  in  inches  depth  over  a  drainage 
area  which  results  from  a  design  rainstorm  is  given  by 

'0: 

Q=        -  V'  (9) 

P  ^  I    +  s 

where  Q  is  accumulated  volume  of  runoff,  P  ~  accumulated  rainfall,  I  is 

a 

an  initial  abstraction,  including  surface  storage,  interception  and  infil- 
tration prior  to  runoff,  and  S  is  potential  maximum  retention  of  water  by 

the  soilo    Each  term  is  measured  in  inches  equivalent  depth  over  the 

watershedc     •    ••   •■    ■  ■  - 


As  a  result  of  studies  of  many  watersheds  in  various  parts  of  the 

United  States „  the  SCS  concluded  that  I  ,  the  initial  abstraction,  can 

a 

be  taken  as  Oo2So    This  causes  equation  9  to  be  reduced  to 


P  +  Oo8  S 


The  potential  maximum  retention  term,        is  usually  replaced  by  an 

expression  involving  a  curve  nimiber,  CNj,  given  by 


CN  .  (11) 

10  +  s 
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The  curve  nianber  CN  is  a  function  of  soil  type  j  land  use ,  vegetative 
covers  soil  temperatures  antecedent  moisture  conditions,  and  other  hydrologic 
soil  characteristics o    With  CN  determined,  equation  9  may  be  solved 
algebraically  J  or  graphically  by  use  of  figure  360  ^ 

The  most  important  factor  affecting  CN  is  the  soil  type^    The  SCS 
has  classified  more  than  ^jOOO  agricultural  soils  in  the  United  States  into 
one  of  h  hydrologic  soil  groups,  A,B,C,  or  Do    Soils  in  group  A  are  highly 
permeable  and  have  low  runoff  potential..    Soils  classified  in  group  D  are 
nearly  impermeable  and  have  high  runoff  potentialo    The  SCS  has  developed 
eactensive  tables  which  incorporate  land  use,  conservation  treatment, 
hydrologic  condition  and  hydrologic  soils  group,  assigning  a  curve  number, 
CNj  to  each  combination »    Antecedent  moisture  conditions  are  accounted 
for  by  arbitrarily  increasing  the  CN  if  conditions  are  abnormally  wet,  and 
lowering  it  if  conditions  are  especially  dry.    Since  most  watersheds  include 
several  different  soil  types  and/or  land  i;^e  practies ,  a  curve  number  is 
established  for  each  subsection  of  the  watershed,  and  a  composite  or  weighted 
CN  is  determined  for  the  basin  as  a  whole  by  prorating  the  areas  of  all 
the  subsectionSc  -       -  -  -. -]  ^  ■   -  ■,. 

Peak  Discharges  from  Basins  Under  2000  Acres 

For  drainage  areas  of  5  to  2000  acres  the  SCS  has  prepared  charts 
showing  peak  discharge  as  a  f motion  of  drainage  area,  curve  number, 
watershed  slopes  geographical  region  and  2i|-hour  rainfall  depth o  Separate 
charts  were  prepared  for  each  curve  number,  for  watershed  land  slopes  of  1%^ 
and  16%^  and  for  two  geographical  regions »    If  a  given  watershed  has  a 
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Figure  -  10. 1  (I  of  2) 


FIGURE  36  —  Soil  Conservation  Service  -  Solution  of  Runoff  Equation 
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different  slope  5,  straight  line  interpolation  is  permitted  using  logarithmic 
graph  papero    Type  1  geographical  region  is  considered  typical  of  maritime 
climates  and  coastal  areas  of  the  United  States  and  Type  2  considered  typical 
of  continental  climates  with  intense  summer  rainstorms o    ESSA  publications 
are  used  to  determine  the  rainfall  depth  corresponding  to  a  2it-hour  duration 
and.  any  desired  storm  frequency <>    Figure  37  shows  the  chart  for  cunre  nmber 
75  9  ^f'  (moderate)  slope,  and  Type  II  stomso  '  r: 

Peak ^Discharge  fran  Gaged  Basins'- 

When  a  basin  under  study  has  streamflow  records  available,  for  which 
unit  hydrographs  can  be  developed ^  the  SCS  uses  these  together  with  r\moff 
volumes  determined  for  design  storms,  in  predicting  peak  discharge  rates  <> 
Base  flow  and  snowmelt  discharge,  if  necessary,  are  added  to  the  discharge 
determined  from  the  unit  hydrographo  /      ■    '  -■ 

Peak  Discharge  f^om  Ungaged  Basins  Using  Normalized  Unit  Hytiro^raphs 

The  SCS  has  developed  112  "normalized"  unit  hydrographs  which  may  be 
used  with  ungaged  watersheds  o    The  112  unit  hydrographs  apply  to  drainage  .  - 
areas  having  times  of  concentration  varying  fVom  lo5  hours  to  72  hours o 
Selection  of  which  unit  hydrograph  to  use  is  based  upon  time  of  concentration 
and  upon  the  ratio  of  accumulative  runoff  at  the  end  of  one  day  divided  by 
total  accumulated  runoff  at  the  end  of  ten  d^s ,  '^2./^]lo°    Time  of  eoneentra- 
tioHj  T^j  for  the  drainage  area  is  determined  as  longest  time  of  bank-full 
channel  flow  plus  overland  flow  timeo    The  direct  runoff  for  one  day  and  for 
ten  days  are  estimated  using  appropriate  curve  number  (CN)  values ^  and  these 
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PEAK  RATES  OF  DISCHARGE  FOR  SMALL  WATERSHEDS 
TYPE  II  STORM  DISTRIBUTION 


SLOPES  -  MODERATE 
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FIGURE  37  — Soil  Conservation  Service  -  Peak  Rates  of  Discharge 
for  Small  Watersheds.     Type  II  Storm  Distribution. 
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estimates  are  reduced  if  the  watershed  area  exceeds  10  square  miles  (areal 
rainfall  is  less  that  point  rainfall),,  and  if  channel  losses  are  significanto 
Channel  losses  are  significant  if  the  stream  is  influento    A  so-called 
"climatic  index,"  Ci,  is  computed  from  the  equation 


100  P 

a 

Ci  s  ^ 


(  12  ) 


a 

where       -  average  annual  precipitation  in  inches 

T    s  average  annual  temperature  in  degress  Fo 

If  the  watershed  has  a  climatic  index  less  than  loOO  the  stream  is  considered 
influent  9  and  channel  losses  are  subtracted  from  direct  runoff  according  to 
tables  prepared  by  the  SCSo    .     j,?   \/u  ■  ■  .■  ,   •  , 

With  the  Q„/Q„^  ratio  and  T    value  established,  the  serial  nimber 
1    lu  c 

of  the  appropriate  unit  hydrograph  is  obtained  from  a  table  o    The  peak 
discharge  is  determined  as  the  maximum  value  of  the  nomalized  unit 
hydrograph  multiplied  by  drainage  area  (in  square  miles)  and  by  Q^^ 
(in  inches  )o    Baseflow  and  steady --state  snowmelt  rates  are  added  if  necessary, 

Peals.    Discharge^  frpm^  Unpaged  Watersheds  Usin^  Synthetie_^Unlt^ _Hydrographs 


A  Soil  Conservation  Service  engineer  j  Victor  Mockus  j  developed  a 
synthetic  imit  hydrograph  method  (reported  in  the  SCS  National  Engineering 
Handbook,  196^4)0    By  Mockus '  method  a  xmit  hydrograph  is  generated  as  the  . 
integral  of  an  "instantaneous  unit  hydrograph,"  the  unit  hydrograph  .  -  „ 

which  theoretically  results  from  a  rainfall  excess  occurring  instantaneously 
on  all  parts  of  a  watershedo    Either  a  curvilinear  or  a  triangular  ■ 
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instantaneous  lunit  hydrograph  may  be  usedo    Except  for  generation  of  the 
synthetic  \mit  l^drograph^  the  peak  discharge  is  determined  in  the  same 
manner  as  in  the  preceding  methods  o 

Project  Formulation^  Prp^ram-Hydrology 

A  priTate  consulting  firm  (CoEdoRo,  InCo )  has  prepared  a  digital 
computer  program  for  the  Soil  ConserTation  Service  which  permits  solu^fcion  of 
the  runoff  prediction  eq,uations  described  above o    The  program  is  versatile 
in  that  it  permits  a  watershed  to  be  divided  into  any  number  of  subareaSj, 
permits  the  'oser  to  supply  a  unit  hydrograph  from  a  gaged  watershed  or 
call  for  the  generation  of  a  synthetic  unit  hydrography  and  pemiits  a  flood 
to  be  routed  through  any  number  of  reserviors  and  any  combination  of 
stream  channels o    A  considerable  amount  of  watershed  data  must  be  supplied 
to  the  computer  before  the  program  can  be  run  for  a  particular  drainage  basin ^ 
but  any  number  of  storms »  having  a  variety  of  antecedent  moisture  conditions  5, 
can  be  solvedo 

The  Soil  Conservation  Service  does  not  usually  run  the  program  with 
actual  storm  datao    Ordinarily 9  hypothetical  design  storm  data  are  lasedo 
In  eases  where  actual  data  are  ti^ed,  first  runs  are  considered  trial.runsg 
and  if  the  computer  output  does  not  yield  runoff  volunes  which  agree  favorable 
with  recorded  volumes  9  adjustments  are  made  in  curve  numbers  or  other  water- 
shed cha?raeteristics  o    Additional  rims  and  adjustments  are  made  until  the 
generated  runoff  volume  satisfactorily  matches  the  recorded  volume o  This 
technique  produces  an  "adjusted"  watersheds  for  which  hypothetical  storms 
can  then  be  run  as  desired,  and  runoff  hydrographs  generatedo 
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TEST  OF  THE  PROJECT  FORMULATION  PROGRAM  -  HYDROLOGY  WITH  DATA  FROM  ONE'''  ' ■ 
STORM  AT  DUCK  CREEK 

A  rainstorm  on  Duck  Creek  watershed  on  June  13-16,  1965  produced  a 
peak  discharge  at  East  Fork  of  205  cfs,  and  at  the  .main  stream  of  Duck  Creek 
of  342  cfs..    The  rain  fell  in  three  principal  bursts,  each  about  24  houi^s 
apart,  with  total  amounts  varying  from  2.26  to  4.30  inches'  in  six  precipit- 
ation stations  on  the  watershed.    Average  total  precipitation  using 
Thiessen  polygon  weighting  was  3.72  inches.    Total  surface  runoff  from  the 
storm  was  0.35  inches.    The  low  yield  is  probably  a  result  of  dry  antecedent 
soil  conditions  and  the  sporadic  nature  of  the  storm.    There  had  been  local- 
ized showers  earlier  in  the  month,  but  not  general-  rain  on  the  watershed  during 
the  three-week  period  perceding  the  June  16  event.  ;  ■  V  /  ■ 

The  storm  was  admittedly  complex,  and  therefore  difficult  to-:modeL., 
After  the  analysis  was  under  way  an  error  in  gage  heights  at  the  Duck  Creek 
gaging  station  was  discovered,  which  meant  that  the  initially  reported  actual 
peak  discharge  and  total  runoff  volume  were  sufficiently  greater  than  the 
true  values  listed  above.  -■  ■  jv,        /  '   ■.  . 

In  order  to  provide  accurate  v/atershed  data  as  input  to  the 
program,  level  lines  were  run  in  the  field  the  full  length  of  each  major  stream 
channel,  with  meander  lengths  and  directions,  numerous  stream  cross  sections 
and  measurements  of  all  stock  ponds  being  obtained  as  well.  ■/^ 

The  first  computer  run  generated  a  runoff  event  having  a  total  yield  of 
0.57  inches.    All  the  curve  numbers  used  were  then  adjusted  downward,  and  a  new 
computer  run  made,  this  time  yielding  a  total  runoff  volume  of  0.43  inches. 
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The  two  runs,  produce^ yhydrographs  having  peak  discharges  of  10P6  cfs 

respectively.  ^  .  tk-a  y^i^r: 

DISCUSSION  OF  RESULTS  •  . 

The  single  storm  on  one  watershed  studied  by  SCS  methods  does  not  permit 
the  drawing  of  any  conclusions  regarding  the  adequacy  of  the  methods  for 
determination  of  peak  discharges.    SCS  hydrologists  would  require  that  the 
program  be  re -run  using  lower  curve  numbers  to  force  the  resulting  runoff 
volume  to  match  that  actually  produced,  before  hypothetical  storms  were  used  to 
produce  design  hydrographs .    It  was  not  possible  to  make  these  additional 
runs,  so  it  can  only  be  speculated  that,  based  on  the  reduction  of  peak 
discharge  from  the  first  to  second  run,  peak  discharge  from  the  "as justed" 
watershed  would  be  reduced  to  820  cfs.    This  is  still  more  than  twice  as  high 
as  the  recorded  peak  flow  of  34-2  cfs. 

The  Soil  Conservation  Service  emphasizes  that  their  methods  are  not 
intended  to  be  used  to  reproduce  actual  runoff  events  that  have  occurred 
(although  the  formulas  and  techniques  used  were  developed  from  actual  storm 
and  runoff  data,  of  course).    The  intention  is  to  produce,  for  a  design 
storm,  the  best  estimate  of  peak  flow  for  spillway  design. 

Interim  Report  #5  (Ferris  and  Williams,  1968)  describing  the  analysis 
of  the  SCS  method  was  reviewed  by  D.  C.  Woo  and  F.  K.  Stovecik,  for  the 
Structures  and  Applied  Mechanics  Division,  Federal  Highway  Administration. 
Woo  and  Stovecik  had  three  comments:    The  SCS  method  is  basically  derived 
for  application  on  agricultural  (as  opposed  to  rangeland)  watersheds;  an 

.aenoni  ■;^il,0  I0  ■■ii"i:Lc'v  ZJ-on'^-i  i.'-'co&  .-.  '^l.'^^J.}^.3^•{l  sm.?:-         ■   ,  f?;:;Bf.:i  jj'/ 
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article  toy  Reieh  and  Hiemstra  (196?)  was  not  eited  in  the  report;  data  from 

the  East  Fork  Dmek  Creek  water  level  recorder  for  the  Jwie  165,  I965  sjtom  should 

be  used  to  extend  the  resialtso    ;  :  -  -.^'V  "  ^  /  .        ■  v-^    .  J: 

The  reTlewers"  first  comment  and  the  results  of  this  investigation  point 
to  the  need  for  extensive  research  to  extend  and  improve  the  SCS  technique 
for  use  on  non-agricultural  watersheds o 

Reich  and  Hiemstra  (I96T)  tested  five  methods  for  estimating  peak  flows, 
using  13^  observed  flood  events  on       streams o    They  pointed  out  that  all-:  -- 
five  methods  tested  are  intended  to  indicate  flood  behavior  which  is  generally 
characteristic  for  typical  conditions  common  in  a  particular  region 3,  rather 
than  to  "predict"  any  single  event o    They  fomid  that  the  SCS  method  under- 
predicted  the  observed  event  88  percent  of  the  timeo     (In  contrast  to  the 
event  studied  in  this  investigation,  for  which  the  method  overpredicted  by 
a  considerable  amount) o    Kent  (1967)9  discussing  Reich  and  Hiemstra's  paper^ 
suggested  that  those  authors  had  used  an  inappropriate  technique  to  compute 
the  estimated  peak  by  the  SCS  methodo    With  estimates  modified  by  Kentg  V": 
the  SCS  method  underpredicted  52  percent  of  the  timeo 

The  East  Fork  Duck  Creek  water  level  records  for  the  June  16^  1965 
storm  were  used  to  make  estimates  of  stream  channel  roughness 9  basin  lag 
time^  etCo    A  separate  computer  run  for  the  East  Fork  alone  could  have  been 
made  by  considering  the  east  fork  drainage  as  a  separate  watershed^  or  the 
East  Fork  hydrograph  could  have  been  used  as  input  to  the  main  watershed 
computer  run  for  the  purpose  of  studying  the  West  Fork  watershed  separately o 
Limitations  on  use  of  the  SCS  computer  precluded  either  of  these  options  being 
triedo  .  ■  ■  . 
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CONCLUSIONS 

The  large  amoimt  of  field  data  that  must  be  obtained  makes  it  appear 
that  the  SOS  methods  are  not  suitable  at  present  for  highway  culvert 
design  except  in  speeial  situations  which  warrant  the  expense  involTed 

in  the  study o 

None  of  the  SCS  methods  assign  return  frequencies  to  the  discharges 
obtainedo    Frequency  is  considered  only  through  the  frequency  of  the  design 
storm  usedg  and  this  may  be  entirely  different  than  the  frequency  of  the  peak 
flow  generated o 

Results  of  this  study  verify  conclusions  which  have  been  reached  by 
others  to  the  effect  that  the  SCS  methods  cannot  be  used  reliably  to  reproduce 
a  hydrograph  from  a  particular  stormc    The  methods  are  intended  to  indicate 
characteristic  flood  behavior  under  typical  conditions  in  the  regiono 

This  chapter  has  reviewed  methods  \2sed  by  the  Soil  Conservation  Service 
for  determination  of  peak  flows,  and  has  described  a  test  of  the  principal 
methods  using  data  from  a  single  rain  storm  which  occurred  on  Duck  Creek 

watersheds 
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CHAPTER  VIT 

MISCELLANEOUS  PEAK  DISCHARGE  METHODS 

VarioTos  public  and  private  organizations  have  developed  peak  discharge 
relationships  for  design  purposes  o    These  are  usually  empirical  formulas 
showing  peak  discharge  as  a  fynetion  of  area,  geographical  factors  and 
sometimes  watershed  characteristics o    Frequency  appears  in  some  of  these 
relationships  as  a  multiplying  factor,  and  does  not  appear  at  all  in  others » 
This  chapter  describes  some  of  these "methods ,  and  shows  a  comparison  of 
results  obtained  by  using  several  methods  with  actual  recorded  peak  discharge 
on  Project  watersheds o    Some  of  the  methods  were  described  in  more  detail  by 
Dodge  (1969)0 

SCS  "COOK"  METHOD  - 

The  Soil  Conservation  Service  published  a  method  in  19^0  for  estimating 
the  i^drograph  peak  flow  Q  frm.  the  equation 

Q  «  P  X  R  X  F  .    "  '  ...  (13) 

for  a  watershed  of  not  more  that  2500  acres  in  sizeo    In  equation  13  P 
is  determined  from  the  chart  in  Figure  380    It  is  a  function  of  drainage 
areaj  and  certain  watershed  characteristics  which  are  described  in  Table  Vo 
^e  rainfall  factor  R  is  really  a  geographical  factor,  and  is  simply 
determined  from  the  map  in  Figure  380    F  is  a  frequency  factor,  taken  to  be 
loOO  for  a  50  year  frequency.    Go 83  for  a  2l+-year  frequency,  and  Go 71  for 
a  10=-year  frequency o    The  rainfall  factor  R  has  not  been  defined  for  Montanao 
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BUREAU  OF  PUBLIC  ROADS  METHOD 

The  Bureau  of  Public  Roads  in  1951  developed  a  method  similar  to  the 

Cook  method  wherein  the  design  peak  discharge  is  computed  from  the  formula 

Q,     o       ^  RF  X  LF  X  LF  FF  X  Q  (ll*) 
design 

Q  is  obtained  as  a  rfunction  of  area  for  basins  up  to  1000  acres  in  size  from 
Figure  39°    A  table  supplies  values  of  LF^  a  land-use  and  slope  factor; 
the  Biap  in  Figure  40  supplies  RFj  a  geographical,  or  rainfall  factor o 
As  with  the  Cook  method  a  frequency  factor  FF  adjusts  the  equation  to 
the  desired  frequency <>    Published  values  of  RF  cover  only  areas  of  the 
United  States  east  of  Montanao 

USBPR  "POTTER"  METHOD 

A  method  for  determining  peak  flows  and  associated  recurrence  intervals 
was  devised  for  the  UoSo  Bureau  of  Public  Roads  by  William  Do  Potter  (l96l)o 
This  method  was  the  result  of  a  research  study  which  as  limited  to  watersheds 
with  areas  less  than  25  square  miles  located  east  of  the  105th  meridian o  The 
study  was  made  on  a  s'ample  of  2ii3  ungaged  watersheds  and  96  gaged  watersheds  o 

The  method  groups  all  watersheds  into  four  physiographic  zones ^  based, 
on  data  supplied  by  the  Soil  Conservation  Service c    These  zones  define  large 
regions  which  are  underlain  by  similar  rock  formations o    The  area  of  Montana 
east  of  the  106th  meridian  is  mostly  located  in  Zone  II ^  with  a  small  section 
located  in  Zone  lo    Figme  \\  shows  a  zone  classification  map  for  the  North 
Central  States o 
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FIGURE  39  ~  Peak  Rates  of  Runoff  for  Drainage  Basins 
Under  1,000  Acres  (from  Bureau  of  Public 
Roads  Manual,  August  1951) 


FIGURE  ^+0  —  Rainfall  Factors.     Use  "with  Figure  39 
in  Estimating  Peak  Rates  of  Runoff 
(from  Bureau  of  Public  Roads  Manual, 
July,  1951) 
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A  topographic  index  T  is  determined  from 


o3L 


(15) 


where  L  -  total  length  of  stream  channel 

and       are  the  slopes  of  the  upper  Oo3  and  lower  OoT  of  the 
stream  channel  respectivelyo 

A  precipitation  index  P  is  determined  to  be  the  amoimt  of  rain  which 
may  "be  expected  to  he  equalled  or  exceeded  on  an  average  once  in  ten  years 

for  a  60  minute  duration o 

A  trial  value  for  the  peak  flow  with  a  10~year  recurrence  interval^ 
^IO(ATP)        obtained  graphicallyo    Figures  h2  and  ^3  are  for  Zones  I  and  II 
respectivelyo    Potter  developed  such  graphs  for  each  of  the  four  zones o 


FigTare  hk  ( for  Zones  I  and  II  or  a  similar  chart  for  zones  III  and  IV) 

A 

is  used  to  determine  an  estimated  value  for  the  topographic  index^  The 
percent  difference  between  estimated  and  measured  values  of  topographic  index 
is  obtained  from 

A 

T      "  T 
AP 


'AP 


X  100 


If  D  is  less  than  30  percent  it  is  assumed  that  the  trial  value  of 

A 

Q^Q  ^j^ipp)  '^^  satisfactory 8  and  this  value  is  accepted  as  the  peak  discharge 
with  10-year  recurrence  intervale    If  D  is  equal  to  or  greater  than  30  percent 9 
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A 

the  ratio  T/T      is  computedg  and  this  value  is  used  in  Figure  k^A  to  obtain 
a  multiplier s  Co    The  trial  value  of  Q^^^  (aTP)       multiplied  hj  C  to  obtain 
the  10~year  discharge o    Potter  suggested  that  peak  discharges  for  other 
recurrence  intervals  could  be  obtained  by  plotting  Q^^  and  Q^^  on  extremal 
probability  paper  and  drawing  a  straight  line  through  the  plotted  points o 

Except  for  the  extreme  eastern  part  of  Montana^  the  Potter  method  is 
not  applicable  in  this  state o    For  purposes  of  comparison  however ^  it  was 
aesumed  that  eastern  Montana  watersheds  could  be  classed  in  Zone  II o  Figure 
k6  shows  the  calculations  used  in  applying  the  method  to  Bacon  Creek  water- 
shedo 

lORTHEM  PACIFIC  RAILWAY  METHOD 

For  many  years  the  Northern  Pacific  Railway  used  the  Ifyers  formula^  wherein 
the  square  feet  of  opening  reqid.red  is  equal  to  a  coefficient  times  the 
square  root  of  the  drainage  areao    Assistant  Chief  Engineer s,  Walter 
Bjorklund  stated  in  a  letter  in  1963  that  they  had  used  coefficients  as  low 
as  Oo5  and  as  high  as  60O0    In  recent  years,  design  has  been  based  upon  records 
of  high  runoff  which  have  been  recorded  in  the  same  general  areao  The 
NF  has  had  several  culverts  which  have  operated  under  a  head,  and  for  which  peak, 
discharge  has  been  computedo    Ustolly these  were  of  sufficient  frequency  that 
they  have  been  able  to  determine  what  to  expect  in  a  particular  areao  By 
superimposing  a  similar  storm  on  other  drainage  areas,    Mro  Bjorklund  indicates 
it  is  possible  to  obtain  the  probable  design  runoffs 
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STEP  A         Identification  of  watershed      BACON  CREEK 


Location:  Topographic  map  quadrangle   Principal  streaxn  MUSSELSHELL 

A-1    Location  of  crossing:    Latitude  ^6  12 '     Longitude    109  50' 
A-2   Area  of  watershed,  A  (planimetered)  =13  .  55 1,000  acres. 

A-3   Principal  stream  length:    L  =  ^^-^  miles;    0.3L  =3-2^  miles;    0.7lJ'^^  miles. 

Elevations  on  principal  stream: 
a,  at  headwateA_965_  ft.;    b,  at  0.7L  above  crossing ^^70  ft^;    c,  at  crossing^^OO  ft^ 

^  (el.  _a  -  el.  b)  -j.  0.3L  =  -k670  )  r  3.2k-  =  195  t3 .  2k  =  60.igt./mi.;  JW^  =  2_.j6 

S   =  (el.  b  -  el.  c)  T  0.7L  =  {kGjO  -k200\  f  7-56  =  ^70  4-7.56  =  62.1^t./mi.;  Vs"  =  7.89 

T  =  (0.3L T      )    (0.7L i  Vs  )  =  (3-24i.7.76)  +  (  7.56^7.89)  =  A 18  +  .957  =  1.375 


STEP  B         From  step  A-1:    Lat.   Long.  

B-1    From  figure  B-la,  b,  c,  or  d:   Zone    II        B-2   From  figure  B-2a,  b,  c,  or  d:  P  =  1 .  0  in. 

STEP  C         From  steps  A  &  B:    Zone     II      A  =13-55    P  =        _     T  =  1-375 

From  figure  C-la,  b,  c,  or  d:    Q  =  1.2    1,000  c.f.s. 

lO(ATP)   


STEP  D 


From  steps  A  &  B:    Zone     jj      A  =13-55    P  =1-0        T  =  1-375 


From  figure  D-la  or  b:    T  =1.200 

AP 


■        Error  =  (T      -  T)  ^  T      x  100  =  ( 1- 20Q  1.  3'^5:  L  20^  100  =  •  ^75^  1.  20^  ik 
AP  AP      ■      •    

STEP_E         From  step  D:    Error  =1^.6 

E-1  &  2         Use  only  if  error  is  equal  to  or  greater  than  30%. 

From  steps  A,  B,  C,  &  D:    Zone  T  =  T      =  Q 

    AP    ^IO(ATP)  

TrT^p  =  f  =   From  figure  E-1:  Coefficient  C  =   

%c)  =  Wp)'"^  =  

E-3         From  steps  C,  D,  &  E-2:    Q    ,         =              Q    ^     =              Error  =  % 
  H      ,    ,  ^IO(ATP)    10(C)     ■ 

When  error  is  less  than  30%   Q     =Q  =  1200  cf  S 

10  ^IO(ATP)   


When  error  is  equal  to  or  greater  than  30%,  Q^^  = 
From  figure  E-3,    Q     =   I8OO  cf  S 


0(C) 


STEP  F         From  step  E-3:     Q     =  ^200     ^     ^  I8OO  cf  S 
  ^  ^10    50   

Plot  Q^^  and  Q^^  on  extremal  probability  paper;  connect  points; 

read  from  extended  curve:     ^     =  I53O    Q      =  2060     Q      =2300     Q      =  268O 

25    -^-^     ^100    200    —   500   

FIGURE  kG.     Potter  Method:     Suggested  Work  Sheet  for  Estimating  Peak 
Rate  of  Runoff  of  Small  Watersheds. 
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Mro  Bjorklund  has  prepared  a  chart  (Figure  ^7)  on  the  basis  of  rmioffs 
which  have  actually  been  experienced  on  the  NPo 

The  line  designated  "Rolling  Wheat  Fields  in  North  DeO^ota  and  Montana 
smd  Rolling  Pasture  in  Minnesota  and  Wisconsin"  was  indicated  by  Mro  Bjorklund 
as  applying  to  much  of  their  Yellowstone  Division  ( f rom  Mandan,  North  Dakota  to 

Livingstons  Montana)  except  not  in  Badlsinds  areas  o 

The  line  designated  "Hilly  Pastures  (Send  Badlands}"  applies  to  most  of 
the  NP  lines  in  eastern  Montana  (Glendive-Sidney „  Glendive-Brockwayg  Glendive- 
Forsythg  Forsyth-Colstrip )  o    Exceptions  are  the  bluff  area  east  of  Miles  Cityj, 
the  bluff  area  east  and  west  of  Rosebud,  and  the  bluff  area  between  Myers  and 
Big  HorBo    In  these  latter  areas  the  line  designated  "Fast  Runoff  (Rugged 
Badlands  and  Rock)"  applies  o 

The  line  designated  "Maximum  Extrene"  indicates  some  extremely  high 
runoffs  which  have  been  encountered  on  lines  other  than  the  NPo    Included  are 
some  floods  in  Missouri  ^  Kansas  and  Nebraska* 

The  W  method  described  is  intended  for  areas  less  that  10^000  acres  in 
Blzm  and  discharges  not  over  IjjOOQ  ef s  ^    Mro  Bjorklmid  indicated^  however,  that 
they  would  be  reluctant  to  design -a-  culvert,  for  a  drainage  area  larger 
that,  3^000  acreas  or  for  a  discharge  of  more  than  1500  efSo 
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UoSo  GEOLOGICAL  SURVEY  METHODS  IN  MONTANA 


Since  inauguratdon  of  the  small-area  peak-flov  highway  program 
in  Montana  in  1955  the  UoSo  Geological  Survey  has  released  three  reports, 
each  presenting  methods  for  obtaining  peak  discharge  rates  having  return 

frequeneies  of  20  or  25  years  o        -  •  -  ,  - 

Berwick  (1958)  proposed  a  method  for  determining  magnitude  and  ■  ■• 
frequency  of  floods ^  applicable  to  any  drainage  area  from  35  to  3,000 
square  miles ^  in  most  of  eastern  Montana,    A  composite  frequency  curve 
expressing  the  ratio  of  the  mean  annual  flood  to  floods  having  return 
period  of  up  to  20  years  was  developed  using  data  from  l6  stations ^  A 
regression  was  derived,  relating  mean  annual  flood  to  drainage  area  and 
mean  elevation  of  the  basino 

An  interim  report  in  I963  by  Boner  revised  Berwick's  report, 
and  expianded  the  relationship  to  cover  drainage  areas  ranging  from 
less  than  1  square  mile  to  3,200  square  miles «    Boner's  technique 
was  the  same  as  that  of  Berwick,  developing  a  composite  frequency 
curve  expressing  ratio  of  mean  annual  flood  to  floods  having  return  ^  . 
period  of  up  to  25  years  using  data  from  19  crest  stage  gages  and 
21  long-term  stations o    A  regression  equation  was  derived,  relating 
mean  annual  flood  to  drainage  area,  mean  basin  elevation,  meander 
length  of  main  stream  course,  and  a  geographical  factorc    The  method  is 
applicable  only  to  areas  of  Montana  east  of  Havre ,  Lewistown  and 
Colimbus,  and  to  basins  under  6OOO  feet  in  elevation » 
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Berwick  and  Boner  each  based  their  procedures  on  a  method 
originally  proposed  in  196O  by  Dalrymple.    Patterson  in  1966  applied 
this  technique  to  most  of  the  Missouri  River  basin  above  Sioux  City,  and 
Bodhaine  and  Thomas  did  the  same  in  196^  for  streams  in  Montana  lying  in 
the  upper  Columbia  River  basin «    Both  of  these  latter  studies  provide 
magnitude  and  frequency  of  expected  floods  for  return  periods  up  to  50 
years  and  for  drainage  basins  larger  than  100  square  miles  o 

Boner  and  Omang,  in  19673  used  the  10-year  flood  as  the  index  flood, 
together  with  ratio  of  25  year  flood  to  10-year  flood,  to  determine  the 
flood  with  25-year  return  frequency o    Montana  was  divided  into  13  geo- 
graphical are^o    The  10-=year  flood  is  a  function  of  which  geographical 
area  the  watershed  is  in 5,  the  area  of'  the  watershed  and,  in  5  of  the  13 
geographical  aresis,  the  average  annual  runoffs    The  ratio  of  25-year 
flood  to  10-year  flood  to  10~year  flood  is  different  for  each  of  the  13 
areas o    Boner  and  Omang  suggest  that  the  reliability  of  their  method 
is  fairj,  and  for  predicting  a  25-year  flood  it  is  pooro 

COMPARISON  OF  VARIOUS  METHODS 

Table  VI  shows  25-year  peak  discharges  obtained  for  the  Project 
watersheds  a  using  the  Boner  s,  Boner-Omang,  Gumbel  and  "Potter"  methods  o 
Design  discharges  for  the  two  watersheds  which  are  under  10,000  acres  were 
also  obtained  using  the  Northern  Pacific  method,  and  these  are  also  shown 
in  Table  VI o    Maximimi  discharges  which  have  actually  been  recorded  on 
the  watersheds  are  shown  for  comparison o 
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The  Boner  .and  the  Boner-Omang  methods  give  25-year  peak  discharges 
which  are  similar  to  the  maximum  recorded  discharges  at  Duck  Creek  and  East 
Fork  Duck  Creek ^  but  they  give  discharges  which  are  smaller  than  those 
which  have  "been  experienced  at  Bacon  Creek,  Hump  Creek  and  Lone  Man  Coulee o 
'The  Gmnbel  method  gives  discharges  which  are  larger  than  any  which  have 
been  experienced  at  Duck  Creek,  East  Fork  Duck  Creek  and  Hump  Creek;  while 
at;  Bacon  Creek  and  Lone  Man  Coulee  the  method  predicts  peak  discharges 
somewhat  smaller  than  those  which  have  been  experiencedo    The  "Potter"  ■ 
method  of  course  was  not  intended  to  be  applied  in  Montana,  and  the  values 
were  included  here  for  comparison  onlyo    The  Northern  Pacific  method  gives 
discharges  which  are  much  larger  than  any  which  have  been  experienced  on 
the  two  watersheds  where  the  method  could  be  appliedo 

True  return  periods  are  not  known  from  the  recorded  discharges  on 
these  streams,  but  it  should  be  noted  that  residents  of  the  Bacon  Creek 
area  report  that  comparable  floods  have  occurred  there  about  once  every  30 
years;  and  that  since  19^8  Lone  Man  Coulee  has  experienced  three  events 
which  were  larger  than  the  25~year  peaks  predicted  by  either  Boner  or 
Boner-Qmango 

COKCLUSIONS 

None  of  the  methods  discussed  in  this  chapter,  appears  to  be  as  reliable 
in  predicting  peak  flows  on  the  five  Project  watersheds  as  the  Gumbel  method 
which  is  a  stochastic  methods 
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CHAPTER  VIII 


COMPARISON  OF  STOCHASTIC  PEAK  FLOW  JffiTHODS 

This  chapter  reports  a  test  of  the  goodness  of  seven  frequency  distributions 
for  prediction  of  peak  discharges  having  specific  return  frequencies  on  the 
five  project  watersheds o     (For  this  test  East  Duck  Creek  was  considered  as  a 
separate  watershed).     The  test  follows  the  procedures  advanced  by  Gupta  (1968) 
in  which  he  tested  ten  methods  on  six  watersheds  located  at  LaCrosse, 
Wisconsin^  ... 

PROCEDURE  ^.  ;  ,^  :  , 

Reference  has  been  made  elsewhere  in  this  report  to  the  multitude  of 
methods  which  have  been  proposed  for  the  prediction  of  the  peak  discharge 
which  has  some  specified  return  frequency »    The  question  which  arises 
is s     Is  one  of  these  methods  more  reliable  than  others  for  use  in  Montana? 
Clearly,  there  is  no  absolute  answer  to  this  question  at  the  present  time 
because  the  exact  correct  discharges  and  frequencies  are  not  known =  However, 
it  is  possible  to  use  statistical  methods  to  compare  various  methods  with 
each  other  and  with  such  recorded  data  as  may  be  available. 

Seven  methods  were  selected  for  comparison,  and  the  recorded  data 
from  the  Project  watersheds  were  used  in  the  analysis.     The  methods  used 
were  the  extreme  value  distribution  (Gumbel's  method) o    The  Log-Pearson  . 
lype  III  method J  four  versions  of  Gringorten's  distributions,  and 
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the  peak  rainfall  frequency-peak  flow  frequency  method  described  in  Chapter 
II  of  this  report o 


The  actual  annual  maximum  discharges  which  have  been  recorded  on  the 
five  watersheds  were  ranked,  and  assumed  return  periods  were  computed 
from  the  widely  accepted  formula 

t    =   (IT 

where  t    -  return  period 

r 

N  ~  number  of  years  of  data 

m  -  ranking  order  of  the  flood  (largest  of  record  -  1,  next  =  2,  etco) 
These  discharges  and  return  periods  were  assumed  for  purposes  of  this  analysis 
to  be  the  absolute  or  correct  values,  so  as  to  have  a  basis  against  which 
to  compare  the  various  prediction  methods » 

For  each  of  the  seven  methods ,  and  for  each  of  the  five  watersheds , 
predicted  discharges  for  the  return  frequencies  established  from  the 
recorded  data  were  determinedo     Each  of  the  methods  is  described  elsewhere 
in  this  report  and  will  not  be  discussed  again  at  this  point  except  to  note 
that  special  calculations  were  necessary  with  the  peak  rainfall  frequency- 
peak  flow  frequency  methodo     Values  of  R^,        and       used  in  this  method  have 
been  established  only  for  i-50  years,  and  accurate  prediction  of  peak  dis- 
charges for  other  return  periods  is  not  possible »    Estimates  of  for 
the  required  return  periods  were  made  for  the  geographic  areas  involved, 
by  using  Gumbel's  technique  to  establish  rainfall  intensities  at  the  lower 
return  periods c  and       were  assumed  not  to  change  at  the  lower  return 
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periods.  Each  of  the  methods  is  purely  stochastic  except  for  the 
peak  rainfall  frequency -peak  flov  frequency  method. 


Results  of  the  study  are  tabulated  in  Appendix  G,  and  are  shown 
graphically  in  Figures  hS  to  52»     (The  peak  rainfall  frequency -peak  flow 
frequency  method  is  identified  on  the  graphs  as  the  ^Robinson'  method.) 


2 

For  each  of  the  seven  methods  the  coefficient  of  determination,  R  ,  was 


then  ccanputed  from 


^2  _  ^         HI  (Predicted  Event  —  Observed  Event)  ^  (l8) 

X3 i Observed  Event  —  Mean  of  Observed  Events) 


The  Gumbel  method  and  all  four  of  the  Gringorten  methods  predicted 
negative  discharges  at  the  lowest  return  f requencies o    For  these  methods 
R"^  values  were  computed  twice,  first  by  using  the  negative  predictions j  > 

and  second  by  assigning  a  discharge  of  zero  to  each  of  the  negative 

2 

predictions o  R    values  are  tabulated  for  the  two  techniques  in  Table  VII » 

2 

The  largest  values  of  R    represent  the  best  correlation  with  the  recorded 


;a;  also  shown  on  Table  VII  are  rankings ,  with  the  method  that 
correlates  best  with  recorded  data  on  a  watershed  being  ranked  l^  the  method 
that  correlates  next  best  being  ranked  2,  etco     Table  VIII  repeats  Table  VII, 
except  that  methods  are  shown  as  columns  and  watersheds  as  rows,  and  the 
rankings  now  indicate  for  which  watershed  each  method  correlated  besto  (For 
example s  the  Gumbel  method  correlated  best  with  recorded  data  on  Duck  Creek, 
next-  best  with  East  Fork  of  Duck  Creek,  etCo  ) 
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FIGURE  hQ  RETURN    FREQUENiCY,  -t^  ,  YEARS  -  ii4  - 


FLOOD  FREQUENCY 

PREDICTED    BY  , 
VARIOUS  METHODS 
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FIGURE  RETURN    FREOUE^'CY^       ,  YEARS         -  115  - 


FLOOD  FREOUENCY 

PREDICTED  BY 
VARIOUS  METHODS 

EAST    FORK    DUCK  CREEK 


—  _^  — — — ■ — — — I  __  

FIGURE  50  RETURKi    FREQUENCr^        ,  YEARS  -  ll6  - 


FLOOD  FREOUEKK 

PREDICTED  BY 
VARIOUS  METHODS 

HUMP  CREEK 


PIGUEE  51  RETURN   FREQUENCY,  t,  ,  YEAR5 


FLOOD  FREOUEN'CY 

PREDICTED  BY 
VARIOUS  MEEHOOi 

LONE  MAN  COULEE 


'  5  10 

FIGURE  52  RETURN    FREQUENCY^  t«  ,  YEARS  -  1^-8 


The  rankings  of  the  seven  methods  on  the  five  watersheds  from  Tables  VII 
and  VIII  were  subject  to  a  Friedman  two-way  classification  test  to  ascertain 
whether  there  are  significant  differences  among  the  methods,  and  whether  there 
are  significant  differences  among  the  response  of  the  watersheds  to  the  various 
fitting  methods  o     The  Friedman  test  involves  the  computation  of  Rank  Chi- 
Sqiiare»  using  the  sum  of  the  ranks  in  each  row,  in  the  formula 

Xr    "    "-TTTT    ?r„^-3r(c+l)  (19) 

rcic+1)  11 

where     "X-T  -Rank  Chi=Square,  to  be  compared  against  tabulated  /L  values 
for  desired  probability  and  degrees  of  freedom  (r-l) 

r  =  number  of  rows 
c  -  number  of  columns 
i\-  sum  of  ranks  in  row  i 

Tables  VII  and  VIII  show  computed  Rank  Chi-Square  values,  together 
with  tabulated  values  for  1%  probability  at  the  specified  degrees  of  freedom. 

DISCUSSION 

Table  VII  shows  that  when  negative  runoff  predicitions  were  assigned 
the  value  zero  the  Gumbel  method  correlated  best  with  recorded  data  on  all 
five  watersheds o    The  Gringorten  methods  generally  ranked  next  hest,  although 
the  Log  Pearson  Type  III  method  ranked  second  best  on  two  watersheds  and  sixth 
best  on  three  watersheds o     The  sum  of  the  rankings  placed  the  Log-Pearson 
method  tied  for  -^th  and  5th  place  with  Gringorten  Type  II  methodo     The  peak 
rainfall  frequency-peak  flow  frequency  method  ranked  poorest,  probably  owing  to 
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inaccurate  evaluat-ion  ofR.,  and  no  adjustment  in  D,  or  F,  for  lower  return 

1  11 

periods  c        '  ' 

When  negative  runoffs  were  permitted  to  remain  in  the  predictions. 
Table  VII  shows  that  the  Gmribel  method  was  hest  on  three  watersheds  and 
hth  best  on  2o     The  Log-Pearson  Type  III  method  was  best  on  two  watersheds 
and  6th  best  on  3=     The  sume  of  the  rankings  placed  the  Gringorten  Type  IV 
method  best,  followed  by  Gumbel,  Gringorten  Type  II  and  Log-Pearson  lype  III 
in  descending  order o     The  peak  rainfall  frequency-peak  flow  frequency  method  was 
again  poorest  o 

All  the  computed  values  for  Rank  Chi -Square  were  larger  than  the 
corresponding  tabulated  Chi-Square  values  at  the  1  percent  probability o  The 
conclusion  therefore  is  that  there  are  indeed  significant  differences  in  the 
correlations  obtained  with  the  seven  methods,  and  also  that  there  are 
significant  differences  in  the  way  the  five  watersheds  respond  to  the  various 
fitting  methods . 

COICLUSIONS 

The  results  reported  herein  indicate  that  of  the  seven  methods,  the 
extreme  value  (Gianbel)  technique  generally  did  the  best  job  of  predicting 
peak  discharges  on  the  five  project  watersheds,  although  when  negative  predictions 
were  permitted  this  method  did  a  poorer  job  of  predicting  flows  at  Bacon 
Creek  and  Lone  Man  Coulee o     The  Log-Pearson  Type  III  method,  which  is  now 
required  for  use  in  frequency  studies  by  all  federal  agencies  did  a  relatively 
poor  job  of  predicting  flows  on  these  watersheds o     The  peak  rainfall 
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frequency -peak  flow  frequency  performed  poorest  of  all  the  methods,  but 
this  is  helieTed  due  to  lack  of  data  from  which  to  adequately  define  the 

terms  in  the  prediction  equation 

In  view  of  the  statistical  difference  found  in  the  five  watersheds  it 
must  be  concluded  that  no  single  method  is  adequate  for  use  throughout 
eastern  Montanao    In  fact,  further  analyses  were  made  using  Duck  Creek  and 
East  Fork  of  Duck  Creek  data.    Since  one  of  these  watersheds  is  a  tributary  of 
the  other  it  might  seem  that  these  watersheds  at  least  would  respond  similarly » 
But  two  tests  5  the  Mann-Whitney  modification  of  the  Wilcoxson  test,  and  the 
Wald-Wolfowitz  Runs  test  both  gave  a  significant  difference  in  the  two 
basins  o 

SUMMARY 

Statistical  tests  to  determine  relative  goodness  of  seven  peak  flow 
prediction  methods  on  five  project  watersheds  have  been  presented  and 

dis cussed o 
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CHAPTER  IX 


STUDY  OF  MUUAL  PEAK  DISCHARGES  FROM  SMALL  DRAINAGE  AREAS 

Annual  peak  discharges  collected  from  228  small  watersheds  by  the  UoSo 
Geological  Sur¥ey  in  cooperation  with  the  Montana  State  Highway  Commission 
ha¥e  been  studied  to  determine  whether  watersheds  in  geographic  proximity 
to  each  other  behave  siirdlarly  hydrologically o    In  those  geographical  areas 
where  such  is  the  case  it  should  be  possible  to  reduce  the  total  number 
of  gaging  stations  required  to  establish  the  hydrologic  properties  of  the 
areao    This  chapter  reports  the  findings  of  this  study c 

BACKGROUND 

The  station-year  method  of  analyzing  rainfall  records  was  proposed  as 
early  as  1936  (see  Linsley,  Kohler^  Paulhus ,  1958) o     The  method  assumes 
that  records  from  several  stations  in  a  limited  area  can  be  combined  and 
treated  as  a  single  record  whose  length  is  equal  to  the  s'um  of  the  individual 
records  o     The  reliability  of  the  analysis  is  determined  by  the  am.ount  of 
dependence  or  independence  between  the  stations  in  the  networks     If  the 
stations  are  so  spaced  that  one  and  only  one  station  measures  each  storm^  the 
data  are  entirely  independent;  but  if  more  than  one  station  measures  each 
stormj,  then  there  are  actually  fewer  independent  stations  and  the  validity  of 
the  method  is  questionable o     The  method  has  been  the  subject  of  considerable 
discussion  and  some  controversy,  but  it  is  generally  believed  to  be  valid 
provided  the  stations  are  truly  independent,  and  located  within  a  meteoro- 
logically homogeneous  areao  ■  ■        ;         ■  1'  . 
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Within  a  hydrologically  homogeneous  area,  and  subject  to  the  same  limit- 
ations on  independence  of  records,  the  method  should  be  equally  valid  when 
applied  to  streamflow  records;  but  no  reference  in  the  literature  was  found 
to  indicate  that  the  procedure  had  been  applied  to  streamflow  datao 

A  recurring  question  in  connection  with  the  UoS.  Geological  Survey' — Montana 
Highway  Commission  cooperative  prograin  to  investigate  the  magnitude  and 
frequency  of  floods  from  sm.all  Montana  watersheds  is  whether  or  not  all  of 
the  more  than  200  gaging  stations  should  be  continued,  now  that  severa.1  years 
of  record  at  each  station  have  been  obtainedo     If  it  could  be  established  that 
two  or  more  independent  watersheds  are  located  in  the  same  hydrologic  unit, 
then  it  should  be  possible  to  pool  the  data  from  all  watersheds  in  the  unit, 
creating  a  record  perhaps  several  times  as  long  as  that  from  any  single  station, 
and  thereby  reduce  significantly  the  error  associated  in  predicting  peak 
flow  magnitudes  and  frequencies o    Within  such  a  hydrologic  unit  it  might  very 
well  be  decided  that  collection  of  further  data  is  not  necessary » 

.     Ao  analysis  of  variance  can  be  made  to  establish  whether  or  not  certain 
watersheds  lie  in  a  single  hydrologic  area,  and  whether  or  not  it  is  valid  to 
combine  data  into  a  single  recordo    This  analysis  cannot,  however,  test  for 
independence  of  recordc 

PROCEDURE 

The  watersheds  included  in  the  cooperative  small  watershed  program  were 
grouped  into  ^8  tentative  hydrologic  units,  using  geographical  proximity, 
same  or  adjacent  drainage  basin ^  and  comparable  elevation  as  criteria  for 
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making  the  groupings o    The  annual  peak  discharges  recorded  from  these  water- 
sheds by  the  UoSo  Geological  Survey  were  inspected  to  ascertain  whether, 
within  each  area,  the  data  were  independent o     The  inspection  showed  that 
frequently  a  single  storm  produced  the  peak  annual  discharge  on  two  or  more 
watersheds o    But  rarely  were  these  discharges  comparable  in  return  period* 
It  is  belie¥ed  therefore  that  with  few  exceptions  the  data  from  the  water 
sheds  can  be  considered  independent  o 

Of  the  228  watersheds  in  the  cooperative  program.,  29  were  excluded  from 
any  group  because  there  were  inadequate  streajnflow  records  available ,  or 
because  they  seemed  to  be  hydrologically  isolatedo    The  hydrologic  groupings 
are  shown  on  a  map  of  Montana,  Figure  53°     Numbers  on  the  map  are  Geological 
Survey  watershed  identification  numbers «    The  watersheds  are  identified  in 
Appendix  Ho  '        '  -.       .  ■    ■  ■  ■  . 

For  each  watersheds  sample  mean  and  standard  deviation  were  computed 
for  the  recorded  annual  peak  flows ;  90^  confidence  limits  were  computed  for  the 
mean  by  using  the  Student's  t  test,  and  an  upper  limit  for  the  standard 
deviation  at  the  90^  probability  level  was  computed  using  a  Chi-Square 
testo     (*These  tests  have  been  described  earlier,  in  Chapter  II} o    Results  are 
tabulated  in  Appendix  Ho 

Within  each  of  the  ^48  geographic  units  shown  in  Figure  53»  recorded 
annual  peak  flows  at  each  station  were  normalized  (reduced  to  flow  per  square 
mile  of  drainage  area)  9  and  the  nomalized  means  and  standard  deviations  were 
subjected  to  an  F  test  to  deteraine  whether  at  the  95  percent  confidence  level 
the  watersheds  are  statistically  different  from  each  other o     (The  F  test  was 
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described  earlierj  in  Chapter  II )o    Results  of  this  analysis  are  tabulated  in 
Appendix  H  and  are  shown  on  Figure  53  where  geographic  units  in  which 
watersheds  are  statistically  different  are  shaded^ 

For  those  areas  which  were  found  statistically  to  be  homogeneous  units, 
pooled  mean  and  standard  deviation  were  computed,  and  these  values  are  tabulated 

for  each  watershed  in  Appendix  Ho 

DISCUSSION 

Twenty  seven  of  the  48  geographic  areas  tested  were  found  not  to  behave  as 
hydrologic  units o     Only  one  of  the  areas  west  of  the  continental  divide  behaved 
as  a  unit;  east  of  the  divide  most  of  the  areas  in  the  Missouri  and  Yellowstone 
headwaters  and  those  in  northern  Montana  behaved  as  hydrologic  \mits  and  data 
from  several  watersheds  can  be  pooled^ 

Those  areas  which  did  not  behave  as  hydrologic  units  might  now  be 
re-examined  dividing  each  area  into  smaller  subdivisions  and  again  testing  for 
hydrologic  similarity;  or  the  areas  might  be  regrouped  placing  some  watersheds 
from  2  or  3  areas  into  one  group  and  testing  for  similarity  <,    A  few  such 
regroupings  were  tried  without  success,  but  a  systematic  or  extensive  program 
of  regrouping  has  not  been  attemptedo 
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FIGURE  53 
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CONCLUSIONS 


Although  this  study  has  not  exhausted  all  possibilities,  hydrologically 
homogeneous  miits  have  been  delineated  in  many  areas  of  the  state c  Within 
these  areas  it  is  belie¥ed  that  recorded  data  are  essentially  independent ,  and 
that  predictions  of  peak  discharge  magnitudes  and  frequencies  may  be  vastly 
improved  by  pooling  data  from  all  the  watersheds.     This  is  evidenced  by  the 
dramaticalj^r  smaller  ranges  on  confidence  limits  for  means  and  the  lowered 
upper  boimd  on  standard  deviation  which  are  shown  in  Appendix  H.. 

It  appears  reasonable  to  conclude  that  to  obtain  the  greatest  benefits 
from  the  cooperative  small  drainage  area  program  in  Montana  it  would  now 
be  desirable  to  relocate  some  of  the  crest  stage  gages  from  watersheds  located 
in  delineated  hydrologic  mits  to  some  of  the  many  areas  of  the  state  where  no 
streamflow  measurements  are  presently  being  obtainedc    At  least  one  or  more 
gaging  stations  within  each  hydrologic  unit  should  be  maintained j  to  obtain  index 
measurements  which  should  be  transferable  to  all  watersheds  in  the  unite  Up 
to  60  gages  might  be  available  for  relocation 

SUMMA.RY 

Twenty-one  hydrologically  homogeneous  units  have  been  delineated 
in  Montana  within  which  it  is  believed  streamflow  data  from  several 
watersheds  may  be  pooled  into  one  record  to  improve  predictions  of  peak 
flow  magnitudes  and  frequencies »    To  maximize  benefits  from  the  cooperative 
small  watershed  program  in  Montana  some  of  the  crest  stage  gages  could  now 
be  relocated  into  areas  of  the  State  where  no  records  are  available o 
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CHAPTER  X 


LENGTH  OF  RECORD  AS  IT  AFTECTS  PEAK  FLOW  PREDICTIONS 

Four  watersheds  which  are  included  in  the  U.S,  Geological  Survey- 
Montana  Highway  Commission  cooperatiYe  small  drainage  area  program  were  studied 

to  see  what  differences  occurred  in  pesk  flow  magnitude  and  frequency 
predictions  as  a  result  of  changing  length  of  recordo     Three  frequency 
distributions  were  used  in  the  study o  :  ^  . 

PROCEDURE 

The  most  frequently  used  peak  flow  prediction  methods  are  usually 
considered  inadequate  when  applied  to  records  of  less  than  10  or  15  years 
in  length 0    To  study  effect  of  record  length  on  predictions,  therefore,  it  was 
necessary  to  select  watersheds  having  streamflow  record  lengths  at  least 
twice  as  long  as  the  acceptable  minimunio     Four  watersheds  which  are  included 
in  the  small  drainage  area  program  were  used  in  this  study.    Each  has  a 
record  length  in  excess  of  30  years o    Three  of  the  four  are  much  larger  in 
area  than  100  square  miles,  but  is  is  believed  that  findings  relative  to  record 
length  on  these  watersheds  will  be  valid  for  smaller  watersheds  as  wello 
Watersheds  selected  are  listed  in  Table  IX o    Watershed  numbers  shown  in  the 
table  refer  to  UoS.  Geological  Survey  identification  numbers;  these  numbers  are 
"used  on  the  location  map.  Figure  53 o 


TABLE  IX o 


Watershed 
Number 

Name 

Area  SqoMi „ 

Record 
Length  yrs  o 

7 

Grasshopper  Creek  near  Dillon 

3h8 

33 

IQ 

McDonald  Creek  at  Winnett 

h21 

31 

138 

Pryor  Creek  near  Billings 

k35 

kk 

209 

Burnt  P'ork  Creek  near  Stevensville 

36 

Three  commonly  used  frequency  distributions  were  used  in  the  study; 
extreme  value  (Gumbel)  method;  log-normal  (Chov)  method;  and  log-Pearson 

Type  III  method o     For  each  watershed  and  each  method,  peak  discharge 
magnitudes  and  frequencies  were  predicted  using  each  half  of  the  record 
separately,,  and  using  the  entire  records    The  results  were  graphed.  Figures  5^ 
to  65 o     Data  and  computations  are  shown  in  Appendix  J\ 

DISCUSSION  .  . 

Examination  of  Figures  5^  to  65  shows  that  no  one  method  gave  superior 
results  on  all  four  watersheds;    and  significant  differences  are  noted 
when  comparing  predictions  made  with  one  half  of  the  data  with  predictions 
made  using  all  the  data^ 

The  Gijubel  method  gave  best  results  of  the  three  methods  tried  on 
McDonald  Creek ^  but  for  the  other  three  watersheds  the  second-half  data  fit 
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the  curve  predicted  from  the  earlier  data  poorly.     On  two  of  the  watersheds 
(Pryor  Creek  and  Grasshopper  Creek)  the  fitting  curve  using  the  first-half 
data  and  the  fitting  curve  using  all  the  data  diverged  greatly  at  the  50-year 
frequency  interval.     On  Burnt  Fork  Creek  the  Gumbel  method  produced  nearly- 
parallel  curves  whether  using  first-half  data  or  all  the  data,  but  on  the  : 
other  three  watersheds  the  corresponding  curves  are  not  parallel,  and  cross 
each  other,, at  approxiinately  the  1.2  to  2.0  year  frequency  interval.  ..  •■: 

The  log-normal  method  was  not  superior  on  any  of  the  four  watersheds.;  ^  -  ■ 
In  each  case  the  second-half  data  fit  the  curve  predicted  from  the  first- 
half  -data; poorly .      The  fitting  curves  using  the  first-half  data  and  all  the 
data  diverged  badl;^''  on  all  watersheds  except  Burnt  Fork  Creek.    With  the 
Exception  of  Burnt  Fork  Creek  the  fitting  curves  crossed  each  other  at     '  >  >'■•■ 
approximately  the  1.0  to  2.5  year  frequency  interval. .i  ^.L  •      j^ri^j        ;  :^">-  ••: 

The  Log-Pearson  Tyve  111  method  was  the  best  method  tried  for  Pryor  Creek; 
however,  on  the  other  three  watersheds  the  second-half  data  fit  the  curve 
from  the  first-half  data  poorly.     On  three  of  the  watersheds,  however,  the 
fitting  curve  ij^ing  the  first-half  data  compared  quite  well  with  the  fitting 
curve  using  all  the  data.     The  curves  diverged  badly  only  on  Grasshopper 
Creek o    On  two  of  the  watersheds,  McDonald  Creek  and  Pryor  Creek,  the  fitting 
curves  crossed  each  other  between  the  3«1  to  3.3  year  frequency  interval. 
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CONCLUSIONS  .    -.rwi  O'  ■ 

Considering  the  fitting  curves  obtained  when  using  the  first-half 
data  as.  compared  to  the  curve  obtained  when  using  all  the  data,  it  appears 
that  the  Log-Pearson  Type  III  method  is  most  satisfactory;  the  divergence 
between  the  curves  at  the  50-year  frequency  interval  was  satisfactory  on 
three  of  the  four  watersheds .     Comparable  curves  using  the  Gumbel  method 
gave  satisfactory  results  on  two  of  the  four  watersheds,  while  the  log-normal 
method  was  satisfactory  on  only  one  watershed. 

Conclusions  from  this  study  are  that  no  single  method  is  satisfactory 
for  prediction  of  peak  flow  magnitudes  and  frequencies  on  all  watersheds; 
and  that  great  discrepancy  can  be  expected  when  predicting  peak  flows  at 
the  50-year  return  period  from  records  as  short  as  10  or  15  years » 

SUMMARY 

Results  of  a  study  of  streamflow  record  length,  using  four  watersheds 
and  three  frequency  distribution  methods  have  been  discussed. 
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CHAPTER  XI 


PHOTOGRMMETRIC  STUDIES 

Reference  has  been  made  earlier  (Chapter  III)  to  the  fact  that  aerial 
photographs  were  taken  of  each  of  the  Project  watersheds  three  times  during 
the  study.     Photos  of  Hump  Creek  and  adjacent  drainages  were  studied  stereo- 
scopically  to  see  whether  such  analyses  might  help  in  evaluating  drainage 
characteristics. 


Several  northward- flowing  tributaries  enter  the  Yellowstone  River  from 
the  south  between  Big  Timber  and  Reed  Point  (about  30  river  miles).  These 
tributaries  mostly  drain  small  watersheds.     Included  are  the  Boulder  River 
(since  it  drains  523  square  miles  it  can  hardly  be  called  a  "small"  water- 
shed). Upper  and  Lower  Deer  Creeks,  Bridger  Creek  (6l  square  miles).  Work 
Creek  (32  square  miles)  and  Hump  Creek  (7.6l  square  miles).     All  these  creeks 
flow  almost  due  north.     East  of  Hump  Creek  there  is  a  sudden  shift  in  direction 
of  the  drainages.    Whistle  Creek  flows  northeasterly,  while  Sectionhouse  Creek 
and  Countrymen  Creek  flow  nearly  east,  as  do  .most  of  the  Yellowstone  tributaries 
from  Reed  Point  to  Columbus,  down  to  and  including  the  Stillwater  River. 
These  drainages  are  shown  on  the  map.  Figure  66,  and  also  appear  in  aerial 
photos  in  Appendix  B  (Figures  B-5  and  B-6). 

The  predominant  geologic  formation  in  the  Hump  Creek  area  is  the  Fort 
Union,  an  extremely  thick  formation  which  contains  beds  of  impermeable  clay 
shale  alternately  with  permeable  sandstone.     The  beds  all  have  a  shallow  dip 
to  the  northeast.     Hump  Creek,  Work  Creek,  and  the  other  streams  to  the  west 
have  breached  one  or  more  of  the  bedding  layers,  and  therefore  have  exposed 
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outcroppings  of  sandstone  and  shale  on  the  eastern  slopes  of  the  water- 
sheds. Headward  erosion  is  occurring  rapidly  on  these  watersheds,  which  are 
quite  steep.  The  streams  to  the  east  of  Hump  Creek  flow  parallel  to  the  dip 
of  the  formation  and  do  not  exhibit  significant  erosion.  Figure  6?  shows  an 
east-west  cross- section  of  the  formations.  Aerial  photos  B-5  and  B-6 
clearly  indicate  the  dip  of  the  formations,  and  show  the  extent  of  the  head- 
ward  erosion.    •-  '■. 

Hump  Creek  and  the  streams  to  the  west  lose  large  quantities  of  water 
to  the  permeable  beds  underlying  them.     Each  watershed  picks  up  small  quan- 
tities of  water  from  springs  on  the  western  slopes  of  the  basins  which  has 
percolated  through  from  the  next  drainage  to  the  west  (Figure  6?)  but  these 
springs  are  small  and  do  not  contribute  significantly  to  surface  runoff. 
Water  which  infiltrates  into  the  permeable  sandstone  layers  is  carried 
under  drainages  lying  to  the  east,  and  probably  under  the  Yellowstone  itself. 
This  subsurface  water  does  not  reappear  as  surface  flow  anywhere  in  the 
vicinity,  and  perhaps  not  even  in  Montana.     Snowmelt  from  the  Absaroka 
mountains  is  captured  by  the  Boulder  River  far  to  the  west.     Water  percolatin 
into  permeable  layers  in  the  Boulder  River  drainage  flows  as  subsurface 
water  far  below  Hump  Creek. 

Streams  lying  east  of  Hump  Creek  i Whistle  Creek,  etc.),  do  not  gain 
significant  quantities  of  water  from  springs,  but  neither  do  they  lose  much 
to  subsurface  flow,  and  hence  they  should  produce  vastly  greater  quantities 
of  surface  runoff  per  square  mile  than  do  Hump  Creek,  Work  Creek,  etc. 
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The  foregoing  discussion  shows  that  geologic  interpretation  by  means  of 
aerial  photographs  provides  extremely  significant  information  regarding 
streatnflow  at  Hump  Creek  and  adjacent  basins.    A  stream  like  Hump  Crek  could 
use  a  much  smaller  culvert,  for  example,  than  Whistle  Creek,  even  though  the 
two  drainages  are  adjacent  to  each  other,  are  nearly  the  same  size,  are  at 
nearly  the  same  elevation,  and  have  numerous  characteristics  of  similarity. 

Similar  analyses  would  seem  to  be  appropriate  in  the  future  as  part  of 
the  hydrologic  investigation  of  any  watershed. 
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•  ,       ;  CHAPTER  XII  '     -     '  ■ 

:.  :  :  CONCLUSIONS  MD  RECOMMDATIONS  :  • ; v-n::;     j c  ;i 

The  research  reported  in  the  preceding  chapters  was  conducted  in  numerous 
phases,  and  the  conclusions  reached  cannot  easily  be  condensed  into  a 
short  space o     The  very  brief  statement  of  conclusions  which  follows  should 
be  supplemented  by  reference  to  the  various  chapters. 

Ic    Precipitation  data  from  Montana  stations  having  long  periods  of 
record  may  be  used  in  a  prediction  formula  to  obtain  estimates  of  peak  dis- 
charges having  a  50-year  return  period.    The  relationship  which  was  derived, 
called  the  peak  rainfall  frequency -peak  flow  frequency  method,  should 
be  refined  by  extension  to  other  return  periods,  and  it  should  be 
checked  by  testing  against  data  from  watersheds  in  other  parts  of  the 
coijintryo     In  its  present  form  use  of  the  method  should  be  limited  to 
furnishing  cross  checks  against  pealc  discharges  predicted  from  other  methods, 

2c    A  multiple -regress ion  equation  which  was  developed  for  five 
Project  watersheds  may  be  used  to  predict  peak  discharge  and  runnoff  volume 
for  certain  rain  or  snow  caused  events «     The  peak  discharge  obtained  has  a 
miknown  return  frequency «     The  equation  has  not  been  tried  on  other  watersheds, 
and  its  reliability  is  not  known »     Future  studies  of  the  hydrology  of  small 
watersheds  for  the  purpose  of  making  peak:  flow  predictions  should  concentrate 
on  the  determination  of  watershed  characteristics  and  collection  of  meteor- 
ological data  as  suggested  in  this  report.     Some  characteristics  and  certain 
types  of  data  seem  to  have  little  bearing  on  peak  discharges  or  rxinoff  volumes » 
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3o     The  Soil  Conservation  Service  method  may  be  used  to  determine  rimoff 
Yolmej  and  to  generate  a  predicted  hydrograph  from  a  hypothesized  rainstorm  or 
from  a  historical  event o    The  runoff  volume  and  hydrograph  obtained 
have  unknown  return  frequencies..    When  historical  data  are  available,  adjust- 
ments may  be  made  in  watershed  characteristics,  thereby  producing  a  runoff 
volume  which  matches  the  recorded  runoff  volume.     The  reliability  of  the 
hydrograph  peak  generated,  however ^  is  doubtful., 

ko     Of  several  multiple  regression  formulas  which  were  tested  (some 
being  essentially  extensions  of  the  ra,tional  runoff  method)  none  were  found 
which  could  be  used  with  much  reliability  to  predict  peak  flows  having  a 
50-year  return  frequency  interval  on  the  Project  Watersheds »     Included  in  this 
group  are  the  Boner,  Boner-Omang,  Potter,  and  Northern  Pacific  formulas. 

5=     Of  several  stochastic  relationships  which  have  been  advanced,  the 
Extreme  Value  (Gumbel)  and  Log-Pearson  T^-pe  III  methods  gave  the  best  results 
when  applied  to  five  Project  watersheds »    Neither  of  these  methods  gave  the 
best  results  consistently o 

60     Certain  areas  of  Montana  have  been  found  to  apparently  be  hydro- 
logically  homogeneous.    Within  these  areas  streamflow  data  from  two  or  more 
watersheds  may  be  pooled  to  produce  a  single  record  which  is  longer  than  that 
from  any  of  the  individual  records e     To  obtain  maximum  benefits  from  the 
current  Small  Drainage  Program^  some  of  the  gaging  stations  which  are  presently 
located  within  these  areas  might  now  be  relocated  in  sections  of  the  state 
where  no  streamflow  records  are  now  being  collected. 


7»     The  length  of  runoff  record  has  a  very  significant  bearing  on  the 
reliability  of  peak,  flows  predicted  by  stochastic  methods  for  high  return 
frequency  intervals o    At  the  50-year  return  frequency  interval,  as  much  as  50 
percent  difference  in  the  predicted  discharge  was  noted  when  using  33  years  of 
record  as  opposed  to  l6c 

Recommendations  are  as  follows:     ;  .  \  •.,   ,  •  y  - 

lo    The  most  reliable  estimates  of  peak  discharges  at  the  50-year  return 
frequency  interval  may  be  made  at  the  present  time  by  the  use  of  stochastic 
methods,  notably  the  Gumbel  and  the  Log-Pearson  Type  III  methods.     If  an   '  " 
estimate  of  the  50-year  peak  flow  is  required  from  a  watershed  where  no  ■ 
streamflow  records  are  available j  the  following  procedure  is  recommended:  " 

a)  Use  the  Gumbel  and/or  Log  Pearson  Type  III  method  on  several 
nearby  gaged  watersheds  and  determine  for  each  the  ratio  of  the  50-year 
peak  discharge  to  the  mean  annual  peak  discharge  or  the  ratio  of  the 
50-year  peak  discharge  to  the  10-year  peak  discharge, 

b)  Use  a  multiple  regression  method  such  as  the  Boner  or  Boner-Omang 
equation  to  obtain  an  estimate  of  the  mean  annual  peak  discharge  or 
the  10-year  peak  discharge  for  the  ungaged  watershedc 

c)  Use  the  average  of  the  ratios  obtained  in  step  a)  as  a  factor  to  be 
multiplied  by  the  mean  annual  peaJ^  discharge  or  the  10-year  peak  discharge 
obtained  in  step  b)^  and  thereby  obtain  an  estimate  of  the  50-year 

peak  discharge  at  the  ungaged  watershed. 
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d)     Use  independent  methods,  such  as  the  peaJs:  rainfall  frequency -peak 
flow  frequency  method  developed  in  this  report j  to  check  the  peak 
flow  obtained  in  step  c)  «     If  major  disagreement  exists,  caution  must 
be  used  in  applying  the  estimated  peak  discharge  to  design  situations. 

2.    Additional  study  of  the  factors  involved  in  the  peak  rainfall 
frequency -peak  flow  frequency  method  is  needed  to  perfect  this  method  and 
improve  its  reliability, 

3o     More  reliable  estimates  of  peak  flow  magnitudes  at  high  return 
frequency  intervals  can  be  made  only  as  longer  historical  records  of  streamflow 
from  small  Montana  watersheds  becom.e  available.    Records  as  long  as  30  years 
must  be  obtained  in  some  areas  of  the  state  before  much  confidence  can  be 
placed  in  estimates  of  50-year  peak  flows.    There  are,  however,  several  areas 
of  the  state  where  the  density  of  gaging  stations  can  now  be  reduced.  The 
gaging  stations  thus  released  should  be  relocated  in  areas  of  the  state 
where  no  records  are  presently  available. 
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